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Abstract 

The F and G form- factors of Ku and the quark condensates are calculated to 0(p 



m 



Chiral Perturbation Theory (CHPT). Full formulas are presented as much as possible. A 
full refit of most of the 0(p 4 ) CHPT parameters is done with a discussion of all inputs and 
underlying assumptions. We discuss the consequences for the vacuum expectation values, 
decay constants, pseudoscalar masses and tt-tt scattering. 
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1 Introduction 

The low-energy realization of QCD is still disconnected from the high-energy structure where 
perturbative QCD applications have made it clear beyond doubt that QCD is the theory of the 
strong interaction. At low energies it is possible to obtain a family of theories within a low-energy 
expansion using only symmetry principles from QCD. This family of low-energy theories is known 
as Chiral Perturbation Theory (CHPT) jlj. The perturbation expansion is in terms of the energy, 
momenta and meson masses, collectively denoted as p. Since this is a perturbative framework, in 
order to compare with experiment it is very important to have calculations performed to as high 
an order as possible. CHPT in the three flavour sector was firmly worked out in Q where all 
parameters were determined to 0(p 4 ). 
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The K e 4 decay 1 , is one of the several processes that can be calculated in the framework of 
CHPT. It has been calculated to tree level, to one-loop Q and to improved one-loop Q. More 
precise K e ^ experiments in progress will go beyond the precision of these calculation. In addition 
K e 4 is the main input in determining three of the CHPT parameters. A lot of work has gone into 
improving the calculation of tt-tt scattering in CHPT, f§, g, g. The precision of the prediction 
on the latter process is now such that these CHPT parameters are the main uncertainty thus 
requiring an extension of the calculation for K e ^ as well. 

This is the purpose of the present work: one further step in the perturbative CHPT series of 
K e 4 or, in other words, to obtain a two-loop result. Using previous two-loop results a complete 
three-flavour and, as far as possible, a model independent result is obtained. We discuss several 
model dependent assumptions, mainly standard, to increase confidence in the generality of the 
numerical results. In addition we also present the calculation of the vacuum expectation values to 
two-loops. 

There are several papers where a complete two-loop calculation is done. Most of them with 
the two-flavour chiral Lagrangian as the effective theory, with the other light flavour and the 
excited states contributing only through the local constants in the CHPT Lagrangian. This two- 
flavour framework seems to be the best one to extract the low-energy parameters of the pion-pion 
scattering, with other states starting to modify the behaviour at higher energies. Also, it is the 
two-flavour sector where chiral symmetry predictions are expected to work best because of the 
clean feature of pseudo-Goldstone bosons for the pions, since the relevant quark mass is so small. 
On the other hand, the two-flavour effective Lagrangian can be applied in a small region of physical 
phenomena. A more interesting Lagrangian has to include the third flavour, but in this case the 
convergence is not so obvious and we can ask the question: is the Kaon mass small enough to 
allow a perturbative theory? To answer this we need a large spectrum of calculations to two-loop 
accuracy. 

At present, besides results to tree and one-loop level, there is some work to try to estimate the 
K e i two-loop contribution Q and there exist several complete two-loop results [ll| for other 
processes. We presented a short summary of the main results with a slightly different input in 
p2[. Here we extend the discussion, present full results and discuss more implications. 



In Sect. 2.1 we present the kinematics with the description of the process and the variables. 



Sect. 2.2 introduces the standard notation for the form-factors followed by a short overview of 



CHPT, Sect. [2^, and a discussion of the available data in Sect. |2.4| . Previous theoretical results 
are described in Sect. [| Next we discuss the calculation, method and checks, Sect. [|and in Sect. 
H some of the 0(p 6 ) constants in CHPT are estimated u sing resonance saturation. Details about 



the 0(p ) constants (LEC), fits and errors are in Sect. I^lj. We also discuss all the input s and 
underlying assumptions there and compare with previous results and a few models, Sect. 3.1.2. 
The resulting theoretical form- factors and partial waves are sho wn in Sect. |6~2]. The comparison 



with a proposed experimental parametrization is done in Sect. 3.2.1. Sect. 6.3 updates the ir-n 



scattering threshold parameters with our fit for the LEC as well as comparisons with available 
low-energy data. Using the LEC fit, we discuss the behaviour of the vacuum expectation values, 



masses and decay constants with the ratio of the quark masses in Sects. 6.4 and 3.5. Finally 
we shortly summarize the results. The appendices provide explicit formulas for the form-factors, 
vacuum expectation values and the one-loop integrals. 



2 Definitions and the Data 
2.1 Kinematics 

We shortly review the variables that parametrize the processes 

K+(p) -> n+ip^TT-ip^e+ipi^M, 

K + ( P ) -> 7T (P+K>-K + (W)^), 

1 The early history of this decay is reviewed in |3|. 
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Figure 1: Kinematical variables in K 14 decays. 

K°(p) -> ir-{p + )^{p^i + {p£)v t {p v ). (2.1) 

In brackets we have labelled the corresponding momentum. For the last two processes we keep 
the momentum notation of the first. The amplitudes for the three processes we denote by , 
T 00 and T~° respectively. Following the original work [|l3| we consider three reference frames, the 
Kaon rest frame, the dipion center-of-mass system [tt — tt plane) and the dilepton center-of-mass 
system (e — v plane). The Ka decays can be parametrized in terms of five kinematical variables 
(see Fig. 0): 

i) Sjr , the squared effective mass of the dipion system. 

ii) S£, the squared effective mass of the dilepton system. 

iii) , the angle between the tt + and the dipion line of flight with respect to the Kaon rest frame. 

iv) 6 1, the angle between the l + and the dilepton line of flight with respect to the Kaon rest frame. 

v) </>, the angle between the tt — tt and e — v planes with respect to the Kaon rest frame. 
These variables can vary inside the range 



Ami < s -k = (P+ +P-) 2 < (wjc - me) 2 , 
mf < s e = (p t + p u ) 2 < (m K - y^) 2 1 

n ^ n n ^ _ n — j. — o_ 



Q<9*,0i<%, O<0<2tt. (2.2) 
Instead of using the previous variables, and besides s OT , cos9 n and 4>, we also use 

U = {p+-p) 2 and u 7r = (p_-p) 2 , (2.3) 



related through 



with 



Stt + t 7T + u 7I = m 2 K + 2 ??4 + s e , 

tn — u n = — 2<7 7r Xcos6' 7r , (2-4) 



X 



-X 1/2 (m 2 K ,s 7T ,S{) , 



A(mi, 7712, TO3) = 77i 2 + 7772 + TT! 2 — 2 (777x7772 + 77117773 + 77727773) . (2-5) 

2.2 Form- Factors 

With the previous notation the amplitude for the decay K + — » tt + tt~ is 

T+- = ^v; s u(pM± - rivbtW - A») , (2.6) 



3 



with 



Vm = -4- 



Ay. = [(p++P-)pF + (p+-p-)pG + (pt+pu)pR]. 

m K 



(2.7) 
(2.8) 



V us is the relevant CKM matrix-element. The other two amplitudes, T~° and T 00 , are defined 
similarly. 

Here we are interested in the F and G form-factors. The H form-factor is already known 
to 0(p 6 ) jl4j and the R form- factor always appears with the mf factor and is negligible for the 
electron case. R is known to 0(p 4 ) Q. The form-factors are functions of s„, sg and cos6* T only, 
or alternatively of s„, t„ and u^. 

The relations between the form-factors and the intensities, easier to obtain from the experiment, 
can be found in |l5| or in || . 

The amplitudes for the three processes of Eq. ( [2.l|) are related using isospin by 



V2 



T 



■00 



(2.9) 



with T 4 - 7 the matrix element defined in Eq. ( |2.6| ). T~~° is anti-symmetric under the interchange 
of the pion momenta while T 00 is symmetric. This also implies relations between the form-factors 
themselves. Observe the different phase convention in the isospin states compared to || where 
M 00 appears with a minus sign because of the Condon- Shortley phase convention. 



2.3 Chiral Perturbation Theory 

CHPT is an expansion in the energy, with the light mesons being pseudo-Goldstone bosons. The 
main feature is the spontaneously broken chiral symmetry. The pseudoscalar fields are the lowest 
mass excitations of the vacuum expectation values. Several realizations are possible, the non-linear 
realization is more useful with the power counting present explicitly. Explicit chiral symmetry 
breaking terms are included as perturbations of the symmetry. For an extensive discussion about 
CHPT we refer the reader to [IJ and references therein. In the following we include the basic 
notation and terms we are working on. 
The pseudoscalar fields are included in 



U(<f>) = u{<}>) 2 = exp( 2 V2$/ J F 1 o) , 



(2.10) 



where 



/ ^ 



V 



V2 V6 



7T 

V2 V6 
K° 



K- 



2r? 8 
V6 J 



(2.11) 



With the chiral symmetry constraints and the transformation properties under this symmetry as 
well as Lorentz invariance and charge conjugation invariance, the Lagrangian starts as 



-effective 
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on 



o-l 



(2.12) 



where the subscripts stand for the power of the low momenta, energies or masses, collectively 
denoted as p. 
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The expressions for the first two terms are (Fq is the pion decay constant in the chiral limit) 

C 2 = ?f{{D^D»U) + (x^U + xt/t)} , (2.13) 

and 

£ 4 = L^D^D^U) 2 + L 2 {D^D V U){D»U^D U U) 

+L 3 (D^D^UD»U^D lJ U) + U{D^D^U) ( X f U + X U f ) 
+L 5 (D^D^U( X i U + U+x)> + Le{x j U + X U f ) 2 
+L 7 ( X ] U - X U ] ) 2 + Lsix^Ux^U + xU^xU^) 
-iL 9 {F^ v D^UD v U^ + F^D^D^U) 

+L w {U^F*UF L n + Hx{F*F R v + F^ v F LfiV ) + H^x) • (2.14) 

The next-to-next-to- leading order £g is a rather long expression and can be found in fl6|| . It 
contains 90+4 terms compared to the 10+2 of £4. 
The notation for the covariant derivative is 

D, L U = d, L U -ir^U + iUl^ (2.15) 

with Ifj, and the left and right external currents related with the vector and axial- vector external 
currents 

r li = v li + a^, l^, = v fi -a fl . (2-16) 

Their field strength tensors are 

F g» =d n»_ff' V i_ i[z",Z"], Fg = 0V -0V -i^.r"]. (2.17) 

Scalar (s) and pseudoscalar (p) external fields are introduced with 

X = 2B Q (s + ip) (2.18) 

where the explicit masses are included through the scalar current s 

s = M + ---. (2.19) 

M. stands for the diagonal quark mass matrix, M. = diag (m u , md, m s ). For later use we define 
the following quantities 

%T = (F£ v + F£ v ) /2 , 
= (F£ v - F£") /2 , 



r„ 



X± = u'x u i U X u > 
V^X = d^X + T^X-XT^. (2.20) 



2.4 Available Experimental Results 

Let us briefly comment on the status of the present data. There are two experiments |l7], on 
which our analysis is based. These two dominate in statistics and precision over previous ones. 
Earlier experiments arc within errors compatible using isospin relations with these two and are 
discussed in 0. In fact our main results rely entirely on Il7j and Q is used to give an idea on 
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the experimental uncertainties involved. In the experiment Jl7| the data are analyzed with the 
partial wave expansion 

F = f s e iS ° + f p e iS "cos6 n + d-wave , 

G = ge iS " + d-wave , (2.21) 

where f s , f p and 3 are assumed to be real. Within the experimental sensitivity no dependence 
of the form-factors on si nor on d-w aves was observed. We confirm this using our calculation as 
shown in Figs. || and [Io| in Sect. 3/2. The expected errors in the planned experiments will allow to 
see the si dependence. The form-factor f p was found to be compatible with zero and consequently 
neglected when the final value of g was derived. As we show below, Fig. [l2|, this is a borderline 
assumption even for the present sensitivity. Neither was there dependence on found for the 
reduced form-factor 

9 = 4-, (2-22) 

Js 

thus the latter is assumed to be constant. 
Furthermore fs{sn) was parametrized as 

fsM = / s (0)(l + A /9 2 ), q 2 = (s v - Ami )/4m*. (2.23) 

Notice that f s is assumed to have only a linear dependence in . As shown in Fig. ^ future exper- 
iments might detect curvature. The linear dependence can only be ensured for s n < 0.11 GeV 2 . 
The assumption of constant g constrains the g form-factor parametrization to be 

g(s„) = g(Q)(l + X g q 2 ) , (2.24) 

with the same slope as the f s form-factor, i.e A = A/ = A g . As is clear, from Fig. |^ the theoretical 
slopes are different but the effect still remain within present errors. Using V us = 0.220 the values 
of ]l7| yield for the form-factors and slope at threshold 

/ g (0) =5.59 ±0.14, 5 (0) = 4.77 ±0.27, A = 0.08 ±0.02. (2.25) 

From |]l8| the same assumptions and using isospin give 

5 (0) = 5.50 ± 0.50, (2.26) 

thus both experiments are marginally compatible. Using the PDG averaging methods |lj| on both 
results we obtain 

3(0) =4.93 ±0.31, (2.27) 

that we use in fit 9. 



3 Form-Factors at Next-to-Leading Order 

For completeness we start with the calculation up to one-loop precision. At tree level -see diagram 
(a) in Fig. |^- both the F and G form-factors are equal and given by a single insertion of £2 [^0) 

F = G =^ 7 (3.1) 

where uik is the physical Kaon mass. At one-loop Q F and G come through the topologies shown 
in Fig. ||, containing two vertices of £2 in (b-c) or a single vertex of £4 in (a). The result can be 
cast in the general form 

F( s „t^) = ^{l + -L|^} + 0( P 6 ), 

G(s w ,t n ,u n ) = ^L-{l + ^E}+0(p 6 ). (3.2) 
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x XK 

(a) (b) (c) 



Figure 2: (a) One-particle irreducible tree level diagram, (b) One-particle irreducible one-loop 
diagrams. Dots refer to strong vertices or current insertions from C 2 , £4 or Cq. External legs 
stand for pseudoscalar or weak current. Internal lines are pseudoscalars only. 



Fx—loop (Gi—ioop) contain three possible contributions: 

i) a unitary correction generated by the loop graphs -see diagram (c) in Fig. |^- given by the 
functions Bi -see App. |Dj- some of which develop an imaginary part responsible for the I = s- 
wave, 1=1 p-wave tt-tt phase-shift. 

ii) A purely logarithmic piece coming from the tad-pole, (b) diagram in Fig. |^. 

iii) A polynomial part which split in two pieces. The first one cures the scale dependence of the 
one-loop functions, making the full process scale-independent while the remaining finite piece is 
tuned to match some experimental observables. 

A straightforward calculation Q leads to 

Fi-ioop = imK-WLl + AL r 2 - 3L r 3 + 8L r 4 + L r 5 - L r 9 ) + 2m 2 K (&L r 2 + 2U 3 - L r 9 ) 
+2«, r (16L; + AL r 3 + LI) - 2U(4L r 2 + 2U 3 - Ug) + 2u„(-\L r 2 + L r 9 ) 
+7/24(~A{ml) + 6A{m 2 K ) + 3A{m 2 )^ + (-3/2ro* + s,)B(m>J, s*) 
+ {ml + m 2 K - t„)/4 B(m\,m 2 K , t„) + {2m\ + m 2 K - u n )/3 B(m 2 , m 2 K , u„) 
/2[B(m 2 K , m 2 K , s w ) + B(m 2 ,m 2 , s T )j + {2m 2 + s„/S) Bi(m 2 , to 2 , s„) 

+ {m 2 + m 2 K /2 - 2U)/3Bi(ml,m 2 K ,U) - u^/'i B Y (m\ , m 2 K , u T ) + s^/2 Bt{m%, m 2 K , s n ) 
+ (-to 2 /3 + 7m 2 K /3 - U)/4Bi(m 2 K , to 2 , t v ) - 2/Ss^ B 21 (m 2 , ml, s v ) 
+ (3m 2 - 3m 2 K - 5/3^)/4 B 21 (ml, m 2 K , t n ) + (3to 2 - 3m| - t v )/4 B 21 (m 2 K , m 2 , t n ) 
-2/3 B 22 {ml , m\ , s*) - 8/3B 22 (ml,m 2 Kl U) - 5/2 B 22 (m 2 K , to 2 , U) , (3.3) 
and for the G form-factor to 

Gi-toop = m*4(-L£ + L r 5 ~ Lg) - 2m 2 K {2U 3 + L r 9 ) + 2s„L r 9 + 2U(AL r 2 + 2L 3 + Ug) 

+2 U7I (-4L r 2 + Lg) + l/24(35A(ml) + UA(m 2 K ) - 3A(to 2 )) 

+ (-ml - m 2 K + t n )/4 B(ml, m 2 K , t n ) + (2ml + m \ ~ u^)/iB{m 2 .,m 2 K ,u n ) 
+ (-ml - m 2 K /2 + 2t 7T )/3B 1 (ml,m 2 Kl t 7r ) - u^/3 5i(to 2 , m 2 Kl u„) 

+ (to 2 /3 - 7/3m 2 K + t„)/4B 1 (m 2 K) to 2 , t») + 3/4(-m£ + m 2 K + 5 /to*) 5 2 i(m 2 , m 2 K , *») 
+ (-3to 2 + 3m 2 . + t^)/4B 21 (m 2 K ,m 2 ,U) - 2B 22 (ml, m 2 , s„) - B 22 {m\ , m%, t w )/3 
-B 22 (m 2 K ,m 2 K , s w ) - B 22 (m 2 K ,m 2 ,ti r )/2. (3.4) 



It is worthwhile to point out at this level some relevant features in Eqs. (3.3), (3.4). For 
instance they depend on the £4 constants L\, L 2 , L 3 , L\, L r 5 and L 9 . Provided we make some 
assumptions (large N c and neglecting d-waves) we can use the low-energy data together with Eqs. 
( |3.3| ), (3.4) to obtain values for L\,L 2 and L 3 . Eqs. ( |3.3| ), (3.4) have also acquired a different 
energy dependence that brings the one-loop result nearer to the experimental values p7| . In fact 
these one-loop corrections are rather large making the low-energy constants evaluation sensitive 
to variations of the treatment of and pseudoscalar masses. 

A possible improvement in the L\ , L 2 and L 3 evaluation is to consider an Omnes-representation 
for both form- factors. In other words, to calculate the dispersive contribution from the next chiral 
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order considering that no intermediate new channels are opened [gj. Even more, renormalization 
group analysis allows to fully calculate the double logarithm contribution . Both showed that 
in this particular case the next order is expected to be quite sizeable. 

This together with the discrepancy between the 577(2) CHPT constants h and 1% obtained from 
the Kn fit || and Roy equation analysis at 0(p A ) ]22| and 0{p & ) |^3| motivated us to calculate 
the two-loop contribution to K^. 



4 Form-Factors at Next-to-Next-to-Leading Order 

4.1 General Technique 

The calculation is done with dimensional regularization (DR) with the extended dimension d = 
4 — 2e. In this scheme the chiral symmetry is preserved and because of the mass independent 
regularization the divergences do not mix different orders in the chiral expansion. 

Once we have all the mathematical tools to calculate all the integrals, the process to obtain 
the form-factors is quite general. We summarize here the main steps and checks. 

One-loop: We reproduce the one-loop result obtained some time ago Q. For reasons explained 
later, in the expansion with the DR parameter e we keep the 0(e 1 ) term. 

Masses and decay constants: A two-loop calculation is a step beyond and with some differ- 
ences with the one-loop calculation. One has to be much more careful in precisely defining 
quantities. E.g. in a one-loop calculation we can safely write all the masses as physical 
because the difference is in the next order. This is the reason why we have to keep in mind 
which are the masses to consider in the two-loop result. The masses can be written as: 

m 2 = m 2 + m 2 ^ +m 2 ^ + ... (4.1) 

with mo the lowest-order mass and where the superscripts (4) and (6) stand for 0(p 4 ) and 
0(p 6 ). For any 0(p 6 ) contribution the masses to,q can be used since the difference is 0(p 8 ). 



However for 0(p 4 ) contribution the masses with the first two terms of Eq. (4.1) should be 
taken. For the tree level the three terms are needed. In CHPT lowest-order and renormalized 
masses are finite and then the result can be in principle be obtained with lowest-order masses 
appearing inside the loops and in the chiral symmetry breaking terms, but we have to use 



the full Eq. (4.1) for the squared momenta. We choose to write all the masses in terms of 
the physical ones. The same discussion can be done for the decay constants; we choose also 
to write the physical decay constant F n instead of the lowest-order one Fq. 

Wave-function- renormalization (WFR): We have to include the WFR for each external leg. 

1/2 1/2 1/2 

That means the matrix element is multiplied by Z^ Z^\ Z V2 , with Zk.tx the renormaliza- 
tion constant for the fields including 0(p 6 ). This is equivalent to include the one-particle- 
reducible diagrams. 

Cancellation of divergences: There are two kind of divergences to cancel. The first one is pure 



polynomial; they cancel with the renormalization of the Ci constants (see Eq. ( 2.12 )). These 
divergences were already found in Q . 

Another kind of divergences are non-local, i.e., can not be expressed as contributions from 
terms in the effective Lagrangian. They appear in the separate contributions and cancel in 
any quantum field theory. But this is a global cancellation, so it taking place, as is the case 
for our calculation, is a strong check of the result. 

A new feature of the two-loop calculations is that the expansion in the DR parameter until 
0(e 1 ) needs to be kept. The reason for that is the existence of terms as LiA(m 2 ) (or as 
LiLj, LiBj, ■ • ■; in general products of one- loop contributions: one-loop x one- loop) where 
A(m 2 ) is a one- loop integral. Once the Li constant is renormalized 

L^Ll + T.X (4.2) 



S 




(d) (e) (f) (g) 



00 (i) (i) 

Figure 3: One-particle irreducible 0(p 6 ) diagrams with only one-loop integrals. The dots are 
0(p 2 ) vertices except the top dot in (d) and (e) which is an 0(p 4 ) vertex. In addition there are 
the diagrams (b) and (c) of Fig. || with one of the dots replaced by a 0(p A ) vertex. 

(h) (i) 

Figure 4: One-particle irreducible 0(p 6 ) diagrams with irreducible two-loop integrals. The dots 
are 0(p 2 ) vertices 

where A = ^- + C and with T and C real and finite numbers, there appears a new finite 
term. 

Double logarithms: Because the nonlocal divergences have to cancel, the double poles in e can 
be obtained from a one-loop calculation. These can be calculated in general plfl or from 
the Li x one-loop graphs for a particular process. We obtain the same result as |2lf for the 
double logarithms and the products Li x Lj and Li x log. 

Three-point one-loop integrals: The contribution from the first four diagrams in the Fig. |, 
(d-g), could be obtained through the renormalization of the meson line in diagrams (b) 
and (c) of Fig. pi Therefore the result should be a product of at most two-point one-loop 
integrals and the naive three-point one-loop integrals from diagrams (e) and (g) must cancel 
when (c) is rewritten in terms of physical masses. 

Irreducible two-loop integrals: The diagrams in Fig. ^ generate sunset-integrals and vertex- 
integrals already known in the literature. In a previous work J9| on the masses and decay 
constants of the pseudoscalar mesons to two-loops, the sunset-integrals are calculated using 
the results obtained in (2^, [2^] . For the vertex- integrals we use another technique explained 
in [^7j and also useful to recalculate the sunset-integrals. This provided another check on 
our integral calculations. Additionally, in the vertex-integrals we obtain the same imaginary 
part from the two-particle discontinuities as using the Cutkosky rules. This check is done 
below the three-particle thresholds. We refer to [£7| for a longer discussion. 

Isospin relations: Isospin symmetry relates the decays as discussed in Sect. [2.2[ . Satisfying 
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these relations is another check of our result. 



Gell-Mann-Okubo (GMO) relation: The last comment in this subsection is related to effects 
from 0(p s ). As explained in the previous paragraph about masses, the correction to the 
0(p 6 ) using the physical masses instead of the bare masses is 0(p s ). Then, to two- loop 
order it is safe to use the GMO relation 3to 2 = 4m|- — to 2 . We can calculate with and 
without this relation. The result is a variation roughly of 10% of the tail of the 0(p 6 ) part 
of the form-factors in Fig. |5|. The plot in the kinematic regime accessible for the experiment 
(s n < 0.16 GeV 2 ) is not modified. This particular effect is mainly due to terms containing 
(3m 2 — Am 2 K + to 2 ) / £„., s„) which appear when reducing the int egra ls to a m ore minimal 



basis. Another example of this type of effects is discussed in Sects. 6.4 and 6.5 



4.2 Results 



Up to 0(p 6 ) we split the contributions as follows 



F^S-r , t-7r , U-K ) — 



m K i Fx-ioop 1 

v^rA + ^r + K i{ LL+ v+ ct 

m K 



V2F V 



I + + i(G u + G V + G ct )} + 0(p s ) 



F 2 



F 4 



(4.3) 



The two first terms have been explained in Sect. ||. The 0(p 6 ) part is split as follows: 
F(G)v: the part from diagrams with vertices from only the lowest-order Lagrangian, it contains 
the pure two-loop unitarity contribution. Fy depends only on the physical masses of the pseu- 
doscalars. 

F(G) c t: is the part of the polynomial contribution depending on the couplings of 

F(G)ll: is the remaining part, in particular all dependence at 0(p 6 ) on the L\ is collected here. 

The diagrams in the Fig. || can be written as product of one-loop integrals and evaluated with 
the standard methods. For the diagrams in the Fig. |^ we profit from previous work with the 
sunset H and the vertex-integrals p7[ ]. The method to manage the sunset- integrals is extensively 
explained in [^| and we refer to this work. For the vertex diagram we refer to [27 since we do not 
explicitly present any formulas containing them. We observe a fairly strong numerical cancellation 
between the contributions from irreducible vertex-integrals and the rest of the terms in Fy and 

Gy. 

We present the results for F(G)ll and F(G) ct in App. A.l and App. A. 2 respectively. For 
F(G)y a significantly longer expression is found. Since it only depends on the physical masses, 
which are not likely to change significantly, a simple parametrization is sufficient. We present it 
for three sets of masses for the three possible decays. Using m#+ and m n + for K + — ► 7r + 7r~ i + vf, 
mx+ and m^o for K + — » 7r°7r°£ + t'£ and to^° and TO^av = ("V+ 

Notice that in all cases we quote F and G as for K + — > tt 
part in cos^^- then gives the form-factors for the other decays. 

The parametrization is given by 



+ rrv0/2 for K° -> n 
~Tr~£ + isi. Taking the even or odd 



F(G) — (<zi + cos0 7r a 2 )(l + sean) + 5^(03 + cos0 w a 4 ) + s 2 (a 5 + cos^a 6 ) 
+s^(a7 + cos 9 n a 8 ) + a^X^ag + cos^aio) • 



(4.4) 



where X = 1/2 ^(m 2 ^ — s v — si) 2 — is^se and a n = yVl — 4m 2 /s w depend on s n and se- The 
values of the parameters are given in Table 0. In Fig. we show how well this parametrization 
fits the calculated expressions. 



5 Estimates of Some G(p 6 ) Constants 

In this section we estimate some of the 0(p 6 ) constants that appear in the results. We assume 
saturation by the lightest vector, axial-vector and scalar mesons, extending the formalism used in 
to the present case. 
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Table 1: The parameters in the fit to Fy and Gy for the three mass cases. 





Fv 


Gv 


F v 


G v 


Fv 


Gv 

v ' V 


TYl OCQpC 


777 t^j. 


777 _l 






771 r,,, 


i n 5 <7i 






1.4238 


0.90367 


1 4^fi 






— il.0993 


+z0.00198 


-i0.88056 


+i0.000468 


-zO.9431 


-zO.000496 


10 6 a 2 


0.1061 


0.060766 


0.09106 


-0.06771 


0.09514 


-0.068565 


10 5 a 3 


-1.5445 


3.0270 


-1.7479 


3.5074 


-1.6368 


3.5262 




+il4.659 


-i0.21235 


+ill.457 


-i0.17035 


+il2.041 


-iO. 15502 


10 7 a 4 


-9.3974 


-5.5732 


-7.8875 


-23.294 


-8.5278 


22.837 


10 5 a 5 


-30.562 


3.4592 


-31.254 


0.037925 


-30.653 


0.035941 




-i0.7821 


+i2.6824 


+il7.946 


+i2.5087 


+il4.810 


+i2.4139 


10 7 a 6 


-0.59204 


2.1912 


-7.1682 


-215.51 


-5.5226 


-206.17 


10 4 a 7 


-3.6231 


-0.803725 


-3.5113 


-0.078793 


-3.6377 


-0.0812 




+i7.6144 


-i0.28955 


+i4.1122 


-i0.26763 


+i4.4861 


-i0.24452 


10 6 a 8 


7.519 


3.1053 


8.3521 


54.23 


8.4695 


50.816 


10% 


-0.54512 


0.69997 


-0.78072 


0.29854 


-0.46175 


0.27801 




+i5.7286 


-i0.02308 


+i6.901 


-i0.29305 


+i6.9589 


-i0.17864 


10Vo 


-7.6007 


-9.7289 


-7.5832 


-9.7041 


-7.6143 


-9.7514 


an 


-2.7610 


-3.0277 


-2.8042 


-3.1412 


-2.7268 


-3.0970 




0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24 

s, s„ 

Figure 5: Fit to the Fy and Gy contribution. Plotted are the fit with lines and the full calculation 
with +. In all cases we plotted cos 6^ = 0,±1 and s e = 0,0.01,0.0225 GeV 2 . 
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With higher orders in the chiral Lagrangian the number of constants to estimate increases. This 
is due to the fact that the chiral Lagrangian is constructed using only some general conditions and 
can describe all the possible theories compatible with these. For any specific theory, QCD in our 
case, there is a large correlation expected between these effective constants. Of course in principle 
they are calculable from the underlying theory but this is not possible at present. The data are 
also not sufficient to determine all of them so we need to estimate using some other method. Their 
main contribution is expected to come from the exchange of higher mass resonances. For these 
we use a Lagrangian with chiral symmetry incorporating extra, higher mass, fields. The precise 
input we use for this is described in this section. 

For the spin-1 mesons we use the realization where the vector contribution to the chiral La- 
grangian starts at 0(p e ). We will only discuss terms relevant for K £4, masses and decay constants. 
Specifically we use for the vector matrix 

" Jv (VK^D+Zx^Kx-D + ^^K,/! (5.1) 



2V2 

and for the axial-vector A 



(5.2) 



with 



and A^ are three-by-three matrices describing the full vector and axial-vector nonets. The rest 
of the notation is as in Sect. [j]. 

There are other terms with the same order in momenta as in the previous Lagrangian, however 
they are related to the anomalous sector, which we do not consider here. 

The scalar mesons are considered through 

Cs = ^(VSV t <S~M 2 s S 2 )+c d (Su^)+c m (Sx + ) + ^(S 2 x+) 

+c 7 (5/+^/D + d r {Sf- la ,fH v ) . (5.3) 

The last two terms are needed for the calculation of the vector and axial- vector two-point functions 
and not interesting for the present discussion. We remark that these are not all the allowed 
contributions from the scalar mesons. We can construct other terms similar to the terms in the 
0(p 4 ) Lagrangian. However that will increases the number of unknown constants to determine 
without improvement for the estimates. It is however one source of error to keep in mind. 

There are two ways to calculate the contribution from these Lagrangians: using the equations of 
motion to integrate out directly the heavy fields; or diagrammatically, considering all the diagrams 
with explicit resonance fields and expanding the result to the 0(p 6 ). In this case it seems faster 
to use the equations of motion. For the spin-1 mesons it is immediate to derive the Lagrangian; 
for the scalars also shifts in their vacuum expectation values have to be considered. 

The result after the integration of resonances is for the vector fields 



^<(V^»)K,X-]>, (5.4, 
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for the axial-vector fields 



f (i) f (i) 

C * = -^(^^-^^^"a) - ^((V (I /^){ U ^ A ^}) , (5.5) 



and for the scalar mesons 



(5.6) 



The three Lagrangians are C(p 6 ). However we remark that the spin-1 Lagrangians are directly 
of this order after the integration. For the scalar case the expansion of the scalar propagator 
produces 0(p i ) as leading contribution but that part is already included via the values of the L\ 
constants. For the work here we take the next term in this expansion to estimate the unknown 
0(p 6 ) constants. 

Besides these ways to estimate the constants it is in some cases possible to obtain the 0(p 6 ) 
couplings from experiment. This has been done in previous calculations |2jJ, where the use of 
finite sum rules could provide the values for the combination of the unknown constants appearing 
in the considered process. It was also done for the pion form-factors [pi] . In all cases the obtained 
value is in reasonable agreement with th e resonance estimate. 



Some of the terms of Eqs. (5^4-5^3) are already directly listed in the full list of the 0(p 6 ) 
Lagrangian terms of fl6}| . The others can be rewritten in terms of them but the resulting list is 
rather long. The combination of terms we are interested in, can be more easily calculated directly 



from Eqs. (5.4-p.q). The resulting resonance contribution is displayed in App. A. 2.1 . Together 



with these we take the couplings 

/y = 0.20, f x = -0.025, 5 y = 0.09, 
a v = -0.014, f A = 0.1, 7a } = 0.1, 
7 ^ 2) = -0.01, c m = 42 MeV, c d = 32 MeV, 

c 7 = 19 ■ 10" 3 GeV -1 , c 7 ~c 7 , (5.7) 

and the masses are the experimental ones flljl . 

my = Trip = 0.77 GeV, m A = m ai = 1.23 GeV, m s = 0.98 GeV. (5.8) 

The values for fy, gv, f x > J a arid c 7 are taken from experiment f28| , |tJ. c m and Cd are from 
comparison with the resonance saturation at 0(p 4 ) pq| . For the remainder we use the ENJL 
model estimates of ||cj . We test the influence of the scalar fields by also performing the fit with 
vector meson contributions only as discussed below (fit 6). 

Before continuing with the next section, it is worth it to summarize some ideas. There are 
90 + 4 0(p 6 ) constants in the most general three-flavour chiral Lagrangian, that are supposed to 
be strongly correlated if they give the low-energy behaviour of QCD. Considering the Lagrangian 
for heavier states, with the same symmetries as the pseudoscalar one, some correlations among the 
constants are obtained. The behaviour of the resonances is different from that for the Goldstone 
bosons: no perturbative theory as for the pseudoscalar mesons exists. However the above procedure 
is arguably the correct picture in the large N c limit. 



5.1 The 7]' Field 

Another possible contribution from heavier states is given through the rj' meson that at low- 
energies does not play as essential a role as for instance the p, except for r\ properties. The rj-rj' 
mixing which affects considerably the r\ properties is known to be fully contained at 0(p 4 ) in the 
low-energy constant Lj. In order to include explicitly the rj' degrees of freedom we can enlarge 
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the pseudoscalar multiplet including the singlet state, in such a way that the r\ and 77' fields are 
a combination of the singlet and octet fields when SU(3) symmetry is broken but SU(2) is still 
respected. An alternative equivalent way is to introduce it as a separate pseudoscalar singlet field 
Pi described by the Lagrangian Q 

£>v> = ld,P 1 d^P 1 l ~M*P* + ^Pi(X-) ■ (5.9) 

At C(p 4 ) the 77' contribution is entirely encoded in L£ and dominates its value. The next order 
expansion in the P\ propagator, 0(p e ), gives 



2 



£>,' = -^-^(x-)cHx-) with d m = 20 MeV. (5.10) 



Similar to the scalar case we can add terms to Eq. (5.9) that will contribute as well to 0(p 6 ). 
The contribution is of relevance to our fit only via the rj mass and the prediction for the 77 decay 
constant, affecting mainly the value for L7 with a small modification for Lg. Again, rewriting Eq. 



(5.10) into the terms of |L6| leads to a long expression so we do not present it. 

We remark that although the two-point functions, used to calculate the decay constants and 
the masses, are obtained with the external currents from the octet, the pole in the propagator of 
course corresponds to the real 77 and the decay constant is the coupling of the real 77 to the octet 
current. 

A more complete list of terms including their integrating out shows only small modifications 
to the phenomenological implications j3lj| . 



6 Phenomenological Applications 

There are two points to address with the two-loop result for Kt&. The first one is to obtain a 
prediction for the form-factors to be compared with experiment and to use our prediction to check 



the validity of assumptions made in the data analysis. This is presented below in Sect. 6.2.1. The 
second one is to extract a fit for the low-energy constants L\ which values will be used to predict 
other quantities of interest. In particular, they are needed to obtain predictions for low-energy 
7r-7r scattering, one phase-shift combination can be obtained independently from the same decay 
data. The comparison with the scattering quantities obtained directly is a consistency check of 
CHPT. We remark that in this two-loop calculation the only way to obtain an imaginary part 
for Kg4 is through the intermediate states of two pions. However the tt-tt amplitude appearing 
here is one-loop. It is then not possible to obtain indirectly a three-flavour two-loop result for 
the rescattering quantities. We therefore use below the known two-flavour two-loop result for the 
scattering || [?]]. 

6.1 Values of the Low-Energy Constants 

The previous values of L\ are obtained from one-loop or one-loop improved calculations [|[ 
These values are obtained through partially independent fits, in particular L\, L^i -^3 are obtained 
from Ku taking fixed the other L\. In the present work we avoid to fix them and perform the 
first global fit of most of the 0(p 4 ) low-energy parameters simultaneously. The total number of 
parameters up to 0(p A ) is 13. They are B m, B m s , F and L\-L\ Q . L r w does not contribute 
to any of the quantities considered here and is thus fully independent. Lg is discussed below. In 
practice we thus need to fit 11 free parameters. The values are updated fitting the results obtained 
for the masses and decay constants to two-loops || and the result of the present work, see also 

The other free quantities in the full result are the C%, that are estimated in the previous 
section, Sect. [|. The values of L\ can be reasonably well predicted from resonance exchange, so a 
similar rough estimate is expected to work for the C\. We test the dependence by removing all the 
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resonance contributions except for the vector ones, whose couplings are determined experimentally 
elsewhere, and the r)', and see if the results change. Changing the values of the C\ by 50% changes 
the fitted values within the fit errors. We also present a reasonable variation of the other inputs. 

Unfortunately, it is not possible to find at present 11 fully independent inputs from direct 
experimental quantities. Below we use as inputs m 2 , m 2 K , m 2 , F^, Fk and three inputs from Kg4, 
/s(Q), g P {0) and A/, defined below. The three additional inputs are the quark-mass ratio m s /rh, 
L\ and Lq. In the future we hope to be able to include additional experimental input to constrain 
these as well. 

The fit is performed by minimizing the \ 2 function defined by 

i=l,6 t =l,6 V ^ X *mput J 

with as weight Aa; m p U ^- the error in the input values £ji n p U t- 

This allows us to fit everything since in all higher orders we have replaced quark masses and 
Fq by the relevant physical quantities. Then we can use the mass and decay constant expressions 
up to 0(p e ) to determine Bqiji and Fq from the physical input afterwards. Let us now discuss the 
various inputs and xt ■ 

Physical masses and decay constants: For these quantities we use the pion masses m,^+ = 
139.56995 MeV, m^o = 0.1349764 MeV, the Kaon masses m K + = 493.677 MeV, m K + = 
497.672 MeV and r] mass m„ = 547.30 MeV, the pion decay constant F n = 92.4 MeV and 
the Kaon decay constant Fk/F* = 1.22 ± 0.01 ||]. 

L4,Lq : The first main feature that we use is the behaviour under the large N c limit. The conse- 
quences for the L\ are discussed in and for the C\ in JTq j . This limit tells us that some 
of the constants are expected to vanish at some unknown scale. This is the case for L\ and 
Lq that are both set zero for the main fit (at m p scale) . Fit 5 in Table |^ assumes them to 
be ze ro at a scale of m 2 instead. For other fits relaxing this assumption we refer to Section 



5.5.1 



lg: Another point is the value for Lg. This constant is related to the charge radius of the pion. 
The value is completely saturated by the p resonance, and for this reason the contribution 
from the logarithms is negligible. Moreover Lg appears with the factor sg in K14, such that 
it is not necessary for the main fit where we use S£ = 0. For the other fits we can safely use 
the value L\ = 6.9 10~ 3 . 

Other If": The remaining constants will be taken as free parameters, at the end six LEC's to fit. 

The quark-mass ratio m s /rh: To 0(p 2 ) it is possible to relate the ratio of the strange quark 
mass m s over the isospin doublet quark mass m = (m u + md)/2, with a combination of the 
meson masses. If we assume that 0(p 4 ) and 0(p 6 ) do not appreciably modify this ratio we 
can use this as input or alternatively the values from QCD sum rules and/or the lattice. 
As standard input we use m s /m = 24 and we check that m s /m — 26 does not lead to 
significantly different results (fit 2 in Table ^) . We remark that there are two possible ways 
to calculate this ratio from the lowest-order meson masses, using or not the r\ mass 



2mL^-mL m a 3 m 2 



"OK 



Ori " l 0ir 



(6.2) 

'Ott 

where mo stands for the lowest-order mass that we calculate in terms of the physical masses, 
and the L\. We included both relations in the fit. The expression of the bare masses in 
terms of the physical masses to two- loop is taken from For the pion mass we use the 
neutral pion mass. For the Kaon mass we subtract the electromagnetic contribution through 
the formula 

m 2 Ksv = \{m 2 K+ +m 2 K0 - 1.8(mJ+ - ™ 2 )) = (494.53 MeV) 2 , (6.3) 
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where the factor 1.8 is a modification due to the corrections to Dashen's theorem |52| 
Whenever we fit or calculate masses we use rn 2 Ka:v and m 2 . The contribution to \ 2 is 



2 



xi+xi=E[ ? ;; • (6 - 4) 

i=l,2 V rh linput / 

Fk/Ftt^ This quantity is calculated using the formulas of Hf], and enters in the fit as 

a ( {Fk/FJ- 1.22 V 

*3 = ^ • ( 6 ' 5 ) 



V o.oi 

IQ4 decay: The expression for the form-factors F and G obtained in the theoretical part is 
compared with the experimental analysis from JltJ . The results we are interested in are the 
extraction of the s-wave contribution from the F form-factor and the p-wave from the G 
form- factor. The others contributing waves are not appreciable ]l7| . This is also confirmed 
in the expressions as discussed below. In |l7j the s-wave / s and p-wave g p are fitted with 
the expressions 

f s (q 2 ) = / s (0)(l + A /g 2 ), 

9 P (q 2 ) = g P (0)(l + \ g q 2 ), (6.6) 

where f s (0) and g p (0) are the form- factors at threshold of the dipion system. The slopes are 
taken to be the same Xf = X g . 

In order to keep the computer time needed in reasonable limits, we use the following ap- 
proximation 

/ s (0)=F(s w ,s^ = 0,cos^ = 0). (6.7) 

We use Sn — (2m 7t + 1 MeV) 2 to avoid numerical problems at threshold. The se dependence 
is very mild, we checked this by performing an alternative fit with ^fsl = 100 MeV (fit 3 in 
Table ||). Only using cos 6^ = corresponds to assume d and higher waves to be negligible. 
Similarly 

g P (0) = G{s v ,s e = 0,cos^ = 0). (6.8) 

We do not make use of the derivative at threshold as the slope. This has both numerical 
problems and does not correspond to the experimental situation. Instead we shift a bit from 
the threshold. The slope is then given by 

_ / |fQ4,^cos^ = 0)| x Ami (f , q] 
} \|F( Sw! ^,cos^ = 0)| 1 ' 

with the value = (2to w + lMeV) 2 and s'^ — (336 MeV) 2 . As an alternative we also 
present a fit with 4 = 293 MeV (fit 4 in Table |). 

These quantities enter in the fit via 

We also present an alternative fit, fit 9, where we include additionally the g p (0) measurement 



in K£ — > tt'tt e + u of Rcf. |18J. Combining errors using the PDG procedure |19| , Eq. (6.1C) 
is modified with 



2 The published version has in App. A. 2 the formulas in terms of lowest-order masses and decay constant, 
contrary to what is stated in the text there. Here we use the formulas in terms of physical masses. 
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Table 2: Results for L^(^l) for the various fits described in the main text. Notice that L£, Lg and 
Lg are input, L\ = 0, -Lg = and Lg = 6.90 10~ 3 . Errors are fitting errors described in the text. 
All L£(/f) values quoted have been brought to the scale /i = 0.77 GeV. The standard values are 
m s /m = 24, ^ = 0.336 GeV and s e = 0. 





Main Fit 


5| 


0(p 4 ) 


fit 2 


fit 3 


fit 4 


fit 5 


fit 6 


fit 7 


fit 8 


fit 9 


10 3 L{ 


0.52±0.23 


0.37±0.23 


0.46 


0.53 


0.50 


0.49 


0.52 


0.45 


0.44 


0.63 


0.65 


10 3 LI 


0.72±0.24 


1.35±0.23 


1.49 


0.73 


0.67 


0.74 


0.80 


0.51 


1.04 


0.73 


0.85 


10 3 LI 


-2.70±0.99 


-3.5 ±0.85 


-3.18 


-2.71 


-2.58 


-2.73 


-2.75 


-2.16 


-2.95 


-2.67 


-3.27 


10 3 U h 


0.65±0.12 


1.4 ±0.5 


1.46 


0.62 


0.65 


0.64 


0.82 


0.67 


1.00 


0.51 


0.60 


10 3 


-0.26±0.15 


-0.4 ±0.2 


-0.49 


-0.20 


-0.26 


-0.26 


-0.30 


-0.26 


-0.23 


-0.25 


-0.26 


10 3 LI 


0.47±0.18 


0.9 ±0.3 


1.08 


0.35 


0.47 


0.47 


0.58 


0.46 


0.52 


0.44 


0.48 


changed 






0(p 4 ) 


m a /rh 


\f$i 




rr.rr 

L14, -kg 


W 


A* 


A* 


5(0) 


quantity 








26 


0.1 


0.293 


-0.3; -0.2 


only 


0.5 


1.0 


4.93 


Unit 










GeV 


GeV 


10" 3 




GeV 


GeV 





With the previous fit procedure we find for the values of the low-energy constants 

10 3 L^ = 0.52 ± 0.23, 10 3 L r 2 = 0.72 ±0.24, 10 s L r 3 = -2.70 ± 0.99 , 

10 3 L£ = 0.65 ±0.12, 10 3 L r 7 = -0.26 ±0.15, l0 3 L r 8 = 0.47 ± 0.18 . (6.12) 

In Table || we present the alternative fits discussed above. The difference with the results of |l^] 
is due to the rf inclusion everywhere, this mainly affects the values of Ly. 

We also included two fits (fit 7 and 8) with a different choice for the scale, n = 0.5 and \i = 1 
GeV. The rescaling of the L\ is considered such that in the table we quote the shifted value at 
fi = 0.77 GeV. These choices imply different scales for the resonance saturation on the C\. As can 
be noticed the changes due to the scale are quite considerable, but almost covered by the fitting 
errors. This suggest that the errors quoted in the main fit are unlikely to be reduced if there is 
more control over the 0(p 6 ) constants. 

The errors are those quoted by MINUIT and are those that give A% 2 = 1 assuming the 
quadratic approximation near the minimum. 

In addition we have performed a fit where V us was changed to 0.226 which showed no significant 
deviations from our main fit. This value is from using unitarity in the CKM matrix and V uc i = 
0.9740 |{n| |§. 

6.1.1 Errors and Correlations 

The problem in determining errors on derived quantities is that the L\ values given above are 
very strongly correlated and that with 6 free parameters the chance that all of them are within 
their one sigma error, is negligibly small. We have thus generated a distribution of sets of L\ that 
have a correlated Gaussian distribution according to the % 2 function defined above. We then keep 
the sets that have a x 2 less than the value corresponding to 68% confidence level. Projecting this 
distribution on the relevant variable leads to larger ranges than the errors quoted above. Those 
are the errors quoted in the remainder since they are more relevant than those calculated from 
propagating the errors of Table |. 

As examples of the correlations we show the plot of L\, L r 2 and L r z together with their projec- 
tions on the two-variable planes and the plot of L7 and L r & in Figure ^]. Each point correspond to 
one set of L\ of the distribution described above. Similar plots can be made for all other combi- 
nations. Notice that this only shows the experimental type of error. Not included are the errors 
due to the restriction to 0(p 6 ) and other theory uncertainties. 

6.1.2 Comparison with Previous Results and Models 

In this section we want to give some insight about the sizeable corrections of the LEC, Eq. ( |6.12 ), 
with respect to the unitary calculation, second column in Table |^, and to show the differences with 
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Figure 6: The sets of L\ within a 68% confidence level range of % 2 . In the first plot we show + 
the values of the first three L\ and x the projections on the coordinate planes. The second plot 
shows the Lj-L r s correlations. 



some existing models based on large N c |34|, [55| . Due to the high correlation between the various 
variables this is not a trivial affair. To make the analysis more transparent we remove first of all 
the constraints on the masses and decay constants, i.e we deal almost with the same assumptions 
as in previous calculations ^ g|. Therefore we fix L r b ,Uj and Lg to their 0(p 4 ) value. All in all 
this amounts to determine L^, L r 2 and L3 from the slope (A), the s-wave and the p-wave (g p ) 
-see Eq. 3J3. The results obtained after that are displayed in fi t A of Table |[ A close look, shows 
that the results are quite similar to those quoted in Eq. ( 6.12| ), thus the constraints imposed by 
the masses and decay constants do not drive at all the differences. To disentangle the pieces from 
the full calculation and compare with || is extremely delicate. For instance, using the dispersive 
method only some pieces of diagrams (h), (j) in Fig. || and of those depicted in Fig. || are kept. 
The rest has no new imaginary part and is therefore modeled by a pure multiplicative constant 
only. Furthermore, the dispersive approach misses the full two-loop renormalization, only the one- 
loop renormalization of F„, masses and WFR is introduced. Besides these obvious remarks a less 
straightforward point needs clarification: in Q the full relativistic resonance propagator is kept, 
thus rcsumming all chiral orders. This has some influence on G. The main point to stress is that 
already at threshold there is a significant shift of the form-factors with respect to the dispersive 
calculation. In Fig. |?] we show this fact, comparing 0(p 4 ), the Omnes improvement and the 0(p 6 ) 
for the F form-factor using the results of Q as input as well as F v = 93.2MeV. Pieces that 
can account for the difference at threshold are the Li x Lj , Li x log, log x log and the resonance 
saturation piece. All of them, except the last one, are 0(p 6 ) contributions not considered in the 
dispersive approach p|. The Omnes improved calculation for F had the right sign but the wrong 
magnitude. In Fig. [7 we have compared the Omnes improved calculation with the full 0(p 6 ). 
Notice the considerable shift already at threshold. To take into account this difference the values 
of the LEC, L\, L\ and L3, have to change substantially. In the figure we have used the values 
of the low-energy constants obtained in |5| , and as a consequence the Omnes calculation matches 
the experimental result, shadow area. With those parameters the 0(p 6 ) result is outside the 
experimental errors. 
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Figure 7: Comparison of the Omnes improved estimate with the full 0(p 6 ) calculation using 
the same old set of values of L\ as input together with F n — 93.2 MeV. The shaded band is the 
experimental result. 



Table 3: Values of L\ for comparison with Eq. ( 6.12] ). The first column is the fit without 
the constraints from the masses and decay constants. Second is the result from the Resonance 
Saturation model. Third and fourth columns are values from other large N c based models as 
described in the text. 





fit A 


Ref. 128 




Ref. | 


34 1 


Ref. 


35, 


38 


10 3 L r 2 
10 3 L r 3 
10 3 L r 5 
10 3 L r s 


0.53 ±0.24 
0.67 ±0.25 
-2.74 ± 1.09 

fixed 

fixed 


0.6 
1.2 
-3.0 
1.4( a ) 
0.9< a > 


0.79 
1.58 
-3.17 
0.43 
0.46 


0.81 
1.62 
-4.24 
1.21 
0.60 



(a) 



Input 
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In our calculation we have assumed a resonance saturation for the 0(p 6 ) coefficients. To some 
extent this assumption was successful at C(p 4 ) p8[ . In the second column of Table || we quote the 
values for this model-approach to some of the LECs. One can see that L\ and Lg in Eq. (6.12) are 
in good agreement with the resonance saturation hypothesis. L 2 is off from this hypothesis, and 
violates the relation 2L\ = L 2 by 20%, which is of the expected size of deviations from the large 
N c limit. As is discussed in ]36j there is some indication that the resonance saturation assumption 
fails in the scalar sector. For this matter, one should compare between different models. It is clear 
that our result is within the range as exemplified by the two examples discussed below. 

We also compare in Table ||| with two model-predictions of CHPT parameters, including also 
In Ref. HU the large N c and chiral limits lead including the leading 

gluonic contribution tc 



4Li = 2L 2 = 



U = — 



96tt 2 



96tt 2 
1 



1 + 0{1/M 6 Q ) 
t: 2 (^GG) 



5N C M% 



0(l/M e Q ) 



(6.13) 



where Mq can be interpreted as a constituent quark mass. Comparing the second column in Table 
| with Eq. fl6.12| ) one can see an improvement in the central value of LJ, L\ and Lg with respect 
to the one-loop result, while the values for L 2 and L^, which essentially agrees with the 0(p A ) 
LECs, went off suggesting the relevance of higher gluonic or large N c corrections. 

The model of 34 was extended to the full ENJL model with gluonic corrections in |j7). The 
resulting relations led to the proposal of a new model, LMD ]3q| . It mimics QCD by taking into 
account the first resonance plus the continuum obtaining the relations 



6Li = 3L 2 



7 3 



4L, 



3 F 2 

8 THp 



15 



8^6 16^ 2 



(6.14) 



where the last equality is obtained using higher orders in QCD sum rules [p8f . Its values are 
given in the third column of Table || and the same comments as in the previous model apply for 
L\,L 2l L^ and Lg, while the central value of L£ turns out to be worse. 



6.2 Kn Form-Factors and Partial Wave Expansions 

We now plot the form- factors F and G at cos 9^ = for the main fit, Eq. ( |6.12 ), and the various 
separate contributions. We show the full result for F in Fig. ||a and for G in Fig. ^b. Notice that 
the convergence of the series is reasonable for both form-factors and that for F the contribution 
from the 0(p 6 ) Lagrangian is fairly small. For G this is dominated by vector exchange with 
experimentally determined parameters so the uncertainty from this source is quite small. 

Let us now look at the dependence on s^. In Fig. || we show the absolute value of F and 
G for ^fsl — 0, 100 MeV and 150 MeV, together with the experimental error band. The lines 
for the higher s^-values stop at the kinematical limits. Notice the change in scale with respect to 
Figures ||a and||b. Given the accuracy of the present data we can safely neglect the dependence 
completely. 

The form-factors F and G can be decomposed in partial wavesQ as |l3|, |l^| 



F(s^,t 7T ,u 7T7 Si) = y^P;(cos6> 7r )/;(s 7r ,s;) - cos 9- IT G{s 7T , t^, U n , Sg) 



i=0 



X 



G(s 7r ,i 7r ,u 7r ,S£) = y^P/(cos si) , 



(6.15) 



3 The predictions for and L% of Ref. |i4[ are sensitive to the precise regularization scheme chosen in the chiral 
quark model. 

4 The partial wave expansion is not simply for F and G since the components with a well denned L z = 0, ±1 
need to be expanded [113, 
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Figure 8: The F (a) and G (b) form- factor at si = and cos^ = 0. The shaded band is the 
result of |l7| . Shown are the full result for the absolute value, and the real part to lowest-order, 
0(p 4 ) and 0(p 6 ). We also show the contribution from the 0(p 6 ) Lagrangian; this is dominated 
by the vector contribution. 
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Figure 10: (a) The partial wave expansion of the combination F = F + a^PL cosd^G / X . Plotted 
are the absolute value of f s and the real part of f s , f p and fd- We also show the f p contribution 
from the second term in F. (b) The partial wave expansion of G. Plotted are the real part of g p 
and gd and the imaginary part of g p . 



where 

P((z) = -^Pi(z). (6.16) 

Pi(x) are the Legendre polynomials and the functions fiis^, se), gi{s^, sg) have as phase the phase- 
shifts 6f (s^) from ir-ir scattering with isospin 1 = for even angular-momentum I, and 1=1 for 
odd I. We perform this expansion numerically for the full F and G expressions. The results are 
plotted in Figs. |(]a and 0b. 

The cZ-waves are very small and the /-waves are not visible on this scale. The contribution 
to f p coming from the term with G, f P (G) is dominant. The difference between f p and f P {G) is 
negligible with the present data as was also found by fiffl. f a has an appreciable imaginary part. 
The imaginary part of f p is also dominated by the part from G. The imaginary part for the other 
waves is negligible. Similarly we see that for G the d- waves are very small. The gf and higher 
waves are not visible on this scale and only g p has a visible imaginary part. 

We can discuss one possible isospin breaking effect. As described for the parametrization in 



Sect. 4.2 we use slightly different masses for the three K e 4 decays. The effect this has on F and 
G is plotted in Fig. [n]. It is obvious that with the present experimental accuracy this source of 
isospin breaking is fully negligible. 

6.2.1 A Parametrization for K + — ► -K + n~e + v 

In |39| two of us proposed a parametrization to be used to extract form-factors and phase-shifts 
from Ka data. Here we check how well the approximations inherent in the parametrization are 
satisfied by our full two-loop calculation. 
The proposed parametrization was 

F = (f s + IX + f" + fiSt)e^ M + f P (T„X cos^e^( s ^ , 

G = (g P + g' v s^ + hsi)e iSlM + ~g d ^X cos 9^ {s ^ . (6.17) 
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Figure 11: The absolute value form- factors calculated with the different masses used for the 
three possible K e 4 decays. The labels indicate the charges of the pion. The curves shown are with 
cosf^ ~ si — 0. (a) |F| (b) \G\. The experimental results are indicated with the grey bands. 

We can now fit these parametrizations to the partial waves obtained in the previous section. 
The tilde notation is because these quantities are related to the partial-wave expansion but not 
identical. 

As foreseen in |3^] the parametrization works extremely well. We have plotted in Fig. [l^ 
how well the absolute value of the s and p-waves of F, f s and f p , and the absolute value of the 
p and d- waves of G are parametrized by Eq. ( |6.17 ). The higher partial waves are completely 



negligible within the experimental precision expected in the near future. The fits are performed 
for Stt < 0.16 GeV 2 only, since that is the region where most events of the data will be. 

6.3 Partial Wave Expansion and Threshold Parameters for tx-tx Scatter- 
ing Amplitudes 

With the main fit result for the 0(p 4 ) low-energy constants, we can discuss how the threshold 
parameters and phase-shifts in 7r-7r scattering are affected with respect to the analysis performed 
in |Q and p3| . For this purpose one expands the ir-ir amplitude as combinations of definite isospin 
amplitudes in the s-channel 

T°{s,t) = 3A(s,t,u) + A(t,u,s)+A(u,s,t), 
T x (s,i) = A(t,u,s) - A(u,s,t), 

T 2 (s,t) = A(t,u,s) +A(u,s,t). (6.18) 
These isospin amplitudes are expanded into partial waves 

oo 

T z (s,t) = 32tt (2/ + l)Pi (cos 6)tf(s), (6.19) 
/=o 

with s,i and u the usual Mandelstam kinematical variables, and 9 is the scattering angle in the 
center of mass frame. Unitary implies that 
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Figure 12: The agreement of the parametrization with our full two-loop result. The points are the 
full result. The lines are the parametrization. We have shown sg = 0,0.01,0.0225 GeV 2 . Notice 
that for comparison purposes, the main parts have been strongly shifted downwards. The fits are 
performed for s v < 0.16 GeV 2 only. 



holds in the elastic region 4to 2 < s < 16m 2 with real phase-shifts 5j (s). s = 4(to 2 + q 2 ) and q is 
the center-of-mass three-momentum of the pions. An expansion near threshold in Eq. ( pT20| ) gives 

1le[t{{ S )]=q 2l {a{ + q 2 b[ + 0{q i )}. (6.21) 

This last expansion defines the threshold parameters, the scattering lengths a\ and the slopes bj . 

The 7r-7r scattering amplitude is known up to 0(p 6 ) -sec also Q for a dispersive calculation- 
in two-flavour CHPT. Expanding this amplitude into partial waves and in q 2 , allows to obtain the 
threshold parameters in terms of the low-energy two- flavour 0(p 4 ) constants, l\. We can study 
the variations of these parameters with the L\ of this work, Eq. ( 6.12| ). However the two-loop 



relations between two- and three-flavour constants is unknown. A full matching would need the 
knowledge of processes involving the relevant constants in both CHPT versions and to the same 
order, this is at present not the case. We use as a first estimate the relations obtained from a 
0(p 4 ) matching 



h 


= 4L[ + 


2L\- 


±v K + <D(p e ), 


'2 


= 4^- 


1 

— VK 

12 


+ 0(p 6 ), 


'.3 


= -8LI 


-ALl 


+ 16L r e + 8L r s ^^ + 0( P 6 


V 


- 8L r 4 + 


ALl- 


\v K + 0( P 6 ) , 



(6.22) 

z, 

with 

H = ^{k + Hf)} aod,p = ^{ln(^) + l}. (6.23) 

We want to emphasize that the relations in Eq. (6.22) are affected by higher -0(p 6 )- chiral 
corrections which might be quite sizeable. Bearing this last issue in mind, in Table || we compare 
previous results and the experimental data (io) with the scattering lengths and slopes obtained 
with the values 

h = 0.40 ± 2.40 , l 2 = 4.90 ±1.0, l 3 = 2.6lt 1 2 9 i , U = 4.05 ± 0.18 . (6.24) 
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Table 4: Threshold parameters in units of 171^+. Set I is l\ = — 1.7 and I2 = 6.1 while Set II is 
Ii = —1.5 and I2 = 4.5. Both sets have 73 = 2.9 and I4 — 4.3. Experimental values are taken from 
pTj| . Errors are from the 68% confidence level ranges only. 





Eq. (3.12) 




0(p 4 ) 


Set I 


C(p B ) Set II 


experiment 




0.219 ±0.005 


0.159 


0.203 


0.222 


0.212 


0.26 ±0.05 




0.279 ±0.011 


0.182 


0.254 


0.282 


0.256 


0.25 ±0.03 


-10a 2 , 


0.420 ±0.010 


0.454 


0.423 


0.420 


0.450 


0.28 ±0.12 


-1062 


0.756 ±0.021 


0.908 


0.755 


0.729 


0.81 


0.82 ±0.08 


10a} 


0.378 ±0.021 


0.303 


0.356 


0.404 


0.38 


0.38 ±0.02 


10 2 6} 


0.59 ±0.12 





0.31 


0.83 


0.56 




10 2 a° 


0.22 ±0.04 





0.16 


0.28 


input 


0.17 ±0.03 


-10 3 60 


0.32 ±0.10 





0.40 


0.24 


0.27 




10 3 a 2 


0.29 ±0.10 





0.32 


0.24 


input 


0.13 ±0.30 


-10 3 6l 


0.36 ±0.9 





0.23 


0.33 


0.22 




10 4 ai 


0.62 ±0.11 





0.20 


0.65 


0.48 




-10 4 6 3 


0.36 ±0.02 





0.15 


0.38 


0.37 





Table 5: The hi quantities as defined in ||. Set I and Set II are defined in the caption of Table 
|]. Errors are the projections of the 68% confidence level domain only. 





Eq. (6.12) 


0(p e ) Set I 


0(p e ) Set II 


10 2 6i 


-6.56 ±2.10 


-9.14 


-8.65 


10 2 6 2 


6.50 ±2.00 


8.88 


8.03 


1O 3 3 


-1.27 ±3.20 


-4.33 


-2.64 


10 3 6 4 


5.44 ± 1.25 


7.08 


4.83 


10 4 6 5 


2.32 ± 1.20 


2.32 


-0.37 


10 4 6 6 


1.16 ±0.23 


1.49 


0.83 



using the Li[ of Eq. ( 6.12 ) and Eq. ( 3.22 ). The errors are given by projections on the relevant 
variable of the 68% confidence level domain. 

In Figure |l^ we plot the predictions of a® and a 2 , for the the set of L\ distributed according to 
the x 2 function as discussed earlier. The projection of the filled ellipse on the axis produces the 
errors in Table ||. 

The interesting result that comes out from Tab le |^ is that with the large difference in the 
input parameters, k, the changes obtained with Eq. fl6.12| ) barely shift the results of set I 0— with 
exception of b\, a® and b®- As was shown in [Q, that is due to the fact that the contribution 
from the double logarithms almost saturates the result to 0(p 6 ). However the new experiments 
are expected to have enough accuracy to distinguish small changes in the scattering lengths and 
slopes. The results of set I and II differ from because of the different value of F n and we 
have kept some 0(p s ) parts which were dropped in the numerics there. Similar remarks apply to 
Figures |l| and [L5|. 

The threshold parameters which were not listed in |Q were calculated by doing the relevant 
expansions in the formulas there to higher order, see App. |^. 

For completeness we also display in Table || the value of the variables bi defined in || as a 
finite combination of the two-flavour low-energy constants and which determines the polynomial 
piece in the n-n scattering at 0(p 6 ). 

In Fig. [l4| we show the difference of phase-shifts 5q — 5{ as a function of the center of mass 
energy. This is the combination measured in decays. As can be appreciated the agreement 
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Figure 13: The prediction for the a[j and scattering lengths for the set of L\ distributed 
according to the correlations of our main fit. Only points within the 68% confidence level domain 
are plotted. 



between Eq. ( p.24| ) and set I is excellent and also with the available data. 

In Fig. [l5| we plot the 1 = 2 s-wave phase-shift Sq as function of the center of m ass e nergy 



together with available experimental data. The phase-shift corresponding to Eq. ( 6.24 ) just 
interpolates between both sets, I and II. 

All in all, the conclusions that can be drawn from the ir-ir scattering analysis is that good agree- 
ment is obtained for all the threshold parameters and with <5q — 8\ and a substantial improvement 
with the data in the 6$ phase-shift case. 

6.4 Vacuum Expectation Values 

The lowest-order effective Lagrangian of QCD is not determined by the pion decay constant 
alone, but involves a second low-energy constant, which is intimately related with the vacuum 
expectation values of the scalar densities in the chiral limit 

(0|gV|0) = -F^B S^. (6.25) 

In the absence of external scalar and pseudoscalar fields it is only observable in combination with 
the quark masses. It is, however, a quantity of theoretical interest and thus deserves study at the 
two-loop level as well. The value of Bq depends on the way the scalar quark bilinear is defined in 
QCD. In particular it depends on the QCD renormalization scale Mqcd- 

The calculation of the one-point scalar Green function without any external pseudoscalar field 
is rather simple as it only involves one- loop x one- loop diagrams at 0(p 6 ). The results are displayed 
in App. |^ in terms of the renormalized quantities, masses and decay constants. The procedure is 



the same as described in Sect. 4.1. In fact, the contributing diagrams form only a small subset of 



those entering in the decay constant evaluation H. This amounts to consider diagrams with one 



vertex involving the external field s part of the matrix \ defined in Sect. 2.3 
We write the full 0(p 6 ) result as 



|W|0) = -F 2 5 {l + <0|gg|0}( 4 > + (0\qq\0) {6) +...,}. (6.26) 
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Figure 14: Difference of phase-shifts 5° — 6\ plotted as a function of the center of mass energy 
-Ejttt. Data points are taken from jL7). The shaded band is the 68% confidence level range. 
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Figure 15: Phase-shift 6q as function of the center of mass energy E^. Data points are taken 
from The shaded band indicates the 68% confidence level range. 
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Note that (0 0) ( ^ 4 - ) and (0\qq\0)^ are scale independent. The soft SU{3) breaking, due to the 
quark masses, introduces an additional dependence on the high-energy constant whose value 
also depends on the precise definition of the quark bilinear in QCD. This ambiguity is forced by 
the need of a subtraction in the scalar two-point function at zero separation. I.e. at the first order 
in perturbation theory 

<0|ot|0) = (0\qq\0) -if d 4 x<0| : (qMq)(x)(qq)(0) : |0) + 0{M 2 ) , (6.27) 

chiral limit J 

which blows up as x — * 0. In contrast [fl2"| the non-analytic contribution is unambiguous. As a 
last point we mention that as QCD predicts Ea| the vacuum condensate increases once the quark 
masses are switched on. 

In order to quote results we need to make certain assumptions on the value of H^. We use the 
scalar dominance assumption that leads to 

HI = 2LI , (6.28) 



with the value of L\ that follows from Eq. (6.12). The numerical results are 



(0\uu\0) = -B F$ [1 + 0.277 + 0.103] , 

(0\ss\0) = -B F^ [1 + 0.832 + 0.247] , (6.29) 

where the numbers correspond to lowest-order, C(p 4 ) and 0(p 6 ) respectively. 

We can also see how the vacuum expectation value changes as a function of the quark masses. 
For this we can use for Bom, Bom s and Fq the values derived using the main fit and the higher 
order formulas expressed in terms of real pseudoscalar masses and F^ . We then insert these values 
into the expressions for the vacuum expectation values in terms of the unrenormalized masses and 



decay constants and vary B^m and Bom s . In Fig. 16 we plotted the va cuum expectation values 
as a function of Tn s /(m s )pY l y S . The difference with the values in Eq. ( 6.29 ) comes from higher 
orders. 



6.5 Masses and Decay Constants 

In this section the pseudoscalar masses are studied using the new values for the LEC, Eq. fl6.12| ). 

All masses are written as in Eq. (4.1) with the higher orders expressed in terms of F^ and the 
physical masses, Jj). The separate contributions for them are 



ml/(ml) phys = 0.746 + 0.007 + 0.247, 
m K/( m K)phys = 0.695 + 0.019 + 0.286, 
m , 2 ,/(^)phys = 0.742 - 0.040 + 0.298. 



(6.30) 



where the numbers are the lowest-order, 0(p 4 ) and 0(p e ) contribution respectively. 

We can now also check how the pseudoscalar masses vary with the quark masses. For this 
purpose we rewrite the pseudoscalar masses as 



in: 



where (m 2 )^ is a function of Bom, Bom s and Fo is defined through 

F^ = Fo(l + F^ +F^ 



(6.31) 



(6.32) 



Plugging in the values following from the lowest-order contribution in Eq. (|6.30 ) we can now vary 
them at a fix ed ratio m s /m — 24. The result is shown in Fig. [l7|a. The difference with the results 
in Eq. ( 6.30 ) is due to higher orders. We believe that the results using the physical decay constant 
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Figure 16: The vacuum expectation values (VEV) as a function of m s / (?ris)phys- Plotted are the 
strange and light quark VEV for the ratio m s /m — 24 and the light quark VEV for fa — 0. The 
curves are normalized to the chiral limit. 




Figure 17: (a) The pseudoscalar masses as a function of the quark mass at fixed ratio fa/ra s . 
The results up to 0(p 2 ), 0(p 4 ) and 0(p 6 ) are shown for the three pseudoscalar masses, (b) The 
pseudoscalar decay constants as a function of the quark mass at fixed ratio fn/m s . The results 
up to 0(p 4 ) and 0(p 6 ) are shown for the three pseudoscalar decay constants. For the r\ it is the 
coupling to the octet axial current. 
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and the physical masses in the higher orders are numerically more reliable since they incorporate 
part of the higher order corrections that we know will appear. 

As has been previously noticed || Eq. ( |S.3C ) the light meson masses at 0(p 6 ) have a 



large contribution which seems in conflict with the convergence of the chiral series. The 0(p i ) 
contribution is much smaller than the 0(p 6 ) one. This is due to two adding effects. First, 
there are cancellations that happen inside the 0(p A ) contribution. E.g. for the Kaon mass the 
contribution is given by 0.066(L^) — 0.047(loops)=0.0f 9 (total p A ). And second and mainly, the 
partial contributions at 0(p 6 ), log x log, Li x log, Li x Lj and pure two-loop, all go in the same 
direction. In other observables we typically had cancellations between these various c ontri butions. 
The corrections to physically measurable quantities are in fact not so large. Eq. ( 6.30| ) can be 
expressed as 



-f = 12.5(p 2 ); 12.7(/); 13.4(/) . 
m 2 



-2 = I6.3(p 2 ); 15.3(p 4 ); 16.4(p 6 ) , (6.33) 
mi 

which shows that the corrections to the ratios are much smaller. 

We have checked that relaxing some of the constraints used in the fits can improve the above 



mentioned behaviour, one of these possibilities is discussed in Sect. |6.5.1| . This does not qualita- 
tively change any of the predictions for decay constants and form-factors made here. Other 
options are ad hoc choosing the values of some of the to diminish the large 0(p 6 ) corrections. 
Some arguments for this can be found in |36| where it is argued that resonance saturation seems 
to fail in the ++ sector. Without new input we feel we cannot discuss further this question. 
For the pseudoscalar decay constants we use the formulas of |j) and obtain using our main fit 

F n /F = 1 + 0.136- 0.076, 
F K /F V = 1 + 0.134 + 0.086, 

F n /F„ = I + 0.202 + 0.069. (6.34) 

The 77-decay constant here is the coupling of the physical rj to the octet axial current. The 
numerical result allows to obtain 

F = 87 MeV . (6.35) 

Using this value of Fq and the lowest-order masses obtained above we can now plot the dependence 
of the decay constants on the quark masses. We therefore rewrite the full formulas in terms of F 



and the quark masses. The results is shown in Fig. |l7|b. Again, the difference with Eq. (6.34) is 
higher order. 

6.5.1 Large N c , Zweig Rule 

In this subsection we relax some of our assumptions to sec whether the unusual relative size of 
the 0(p 6 ) versus the 0{p A ) contribution in Eq. ( 6.30| ) can be avoided. 



We have performed fits with m s /m = 20 and m s /m = 30 but otherwise the same input as the 
main fit. Neither of these changes the relative size of the contributions to the masses qualitatively. 

Another possibility is to relax the accepted common wisdom about the large N c (44| behaviour 
of the low-energy constants. In particular we will focus on L\ and Lg, recapitulating some of the 
reasoning of g] . 

In principle the main question is whether the corrections of the relative order m s and to 2 to 
the ground state are small in comparison with the scale of the theory. If the answer is no, clearly 
one can not trust large N c arguments in CHPT. Large corrections would also cast doubt on CHPT 
beyond one-loop calculations in the three-flavour sector. 

At 0(p 4 ) the effects of m s on (0|uu|0), F T and m 2 are small if both L\ and Lg ~ 0. The main 
question is now if higher chiral orders in those quantities are still protected against a strong m s 
dependence. As it is shown, there is a range in L\ and Lg values that allows large N c and CHPT 
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at 0(p 6 ) to coexist without conflict. Thus there is no reason to argue for large deviations from 
the Zweig rule. 



Setting U = 1 in Eq. (2.14) the large iV c counting comes automatically, after using some trace 
identities || 

0{N C ) Li, L 2 , L 3 , L 5l L$, L g , L Wl Hi and H 2 , 

0(1) 2L 1 - L 2 ,L 4 , and L 6 . 



The behaviour of the constant L-j depends on the treatment of the singlet rji (|, 45 . In the con- 
stants of 0(1) there are two different behaviours: i) the combination 2L\ — L 2 is scale independent, 
but ii) L\ and Lg are scale dependent. This last dependence should cancel via the meson loops 
that explicitly violate large N c . In order to estimate possible sizes we can just vary the scale inside 
a reasonable range. Using 



IALI 



\m^ 2 )-L r M\ 



167T 2 



■hi 



with the choice [i\ '. [i 2 — 2, leads to 



(0 ± 0.5) • 1CT 3 and L r 6 = (0 ± 0.3) • 10" 



(6.36) 



(6.37) 



The central value correspond to the main fit presented in Table ||, i.e. considering that large N c 
is reached at t he p scale. Fit 5 corresponds to use large N c at the rj scale. Neither Eq. ( 6.12 ), as 
shown in Eq. ( |6.30 ), nor fit 5 have the expected ordering of 0(p 4 ) and 0(p 6 ) contribution to the 
pseudoscalar masses. 

Other constraints on L\ and -Lg are from a recent analysis [36| based on the positivity of the 
fermionic measure giving a lower bound 10 3 Lg(M p ) > 0.13 and on QCD sum-rules jl6| which 
allow a band 0.2 < 10 3 L r 6 (M p ) < 0.6. The latter reference also quotes W 3 Ll(M p ) ~ 0.4 with 
unknown errors. The values and bounds from [ |36"[ |46] | use 0(p 4 ) expressions so these ranges might 
well change when those calculations are performed to 0(p 6 ). 

Let us now check for which ranges of L\ and L§ we find results where all the quantities 
considered up to now have an expansion of the form 



X = X 2 +X 4 +X 6 with \X 2 \ > \X 4 \ > \X 6 



(6.38) 



where the subscript refers to the chiral order. This is for reasonable variations always the case 
for f s (0), ffp(0), A/, (0|ss|0) and F v /F n . The ratio s of t he pseudoscalar masses a lway s have small 
corrections even if they do not always obey Eq. ( 3.38 ). We will refer to Eq. ( |6.38| ) above as a 
convergent series. For (0|um|0) H 2 has only a small influence but it is very large for (0|ss|0). For 
this quantity we thus allow \X 2 \ < \X 4 \ since this can easily be changed by varying H^. 

The problematic quantities are F^/Fq, Fx/F^, the pseudoscalar masses separately and the 
vacuum expectation value (OlwulO). 

In Fig. [l^ we show the allowed region in the plane L\ — L r & consistent with Eq. ( 6.38 ) for all 
quantities^]. The allowed bottom right region extends further to the right but with ( Fk/F ^a < 
(Fk/F^q even though both are small. As can be seen there is a region where Eq. (6.38), large 
TVc constraints and the other bounds are reasonably satisfied, taking into account possible higher 
order corrections to the bounds. Notice that a large deviation from the Zweig rule is not needed 
thus the bounds, Eq. ( |6.37| ), still hold. 

Let us mention some general features in the behaviour of these quantities as a function of L\ and 
ig. For instance, the decay constant has the broadest range (in both parameters) of convergence, 
while the mass seems to be the most restrictive. A low, negative value for both L\ and Lg is 
penalized for the mass and the vacuum expectation value. If we increase L\ fixing the value of Lq 
the vacuum expectation value recover once more the convergent pattern. The mass behaviour is 
quite unexpected: large, absolute values of Lg reorder the series making 0(p 4 ) > 0(p G ), but then 
0(p 4 ) starts to compete with 0(p 2 ) contribution. 

5 Quantities are only calculated on the points of the grid. 
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Figure 18: The L\-U§ plane. The grey circles arc the points where Eq. (6.38) is satisfied for 
all quantities considered in this manuscript. The grey band shows the bounds on Lg from QCD 



sum-rules. The dash-dotted line is the lower bound that follows from positivity constraints |3q| 
and the thick dashed line the value of L\ of Ref . . Fits were only performed for points on the 
grid. 



7 Summary 

Let us summarize our results. We have described the calculation to 0(p 6 ) in CHPT of K14 decays 
and of the vacuum expectation values, giving some general lines on the methods that we used 
and all the various checks on the final results. The explicit formulas have been presented as 
much as possible. The long expressions involving the vertex- and sunset-type integrals we simply 



parametrized, Eq. (4.4). The rest of the expressions, involving the low-energy constants, are fully 
listed in App. [A|. 

We have assumed resonance saturation to estimate the new 0{p 6 ) parameters that appear. 
With all this at hand, and large N c assumptions on L\ and Lg we updated most of the values of the 



0(p 4 ) parameters, L\, Eq. Q6.12 ), testing the validity of the large N c relation (U 2 — 2L^)/Lg ~ 0. 



To this end we made use of the existing K14 data, pseudoscalar masses and decay constants. We 
varied the inputs over a fairly wide regime and discussed in detail the underlying assumptions 
and their consequences once some of them are relaxed, see Table |J Among the assump tions it is 



worth to mention the discussion on relaxing the large N c behaviour of L\ and Lg, Sect. p.5.1| . We 
show the allowed region in the L^-Lg-plane, where CHPT, as a well defined perturbative series, 
can exist without conflict with any of the known quantities up to 0(p 6 ). 

We also present the comparison with previous results and model estimates, pointing out the 
source of the difference. As it becomes immediately clear the main source of the theoretical 
uncertainties resides in the 0(p 6 ) constants. We mention briefly how to treat the errors and show 
one way of treating the correlations using a 68% confidence level distribution of the L\. 

Once the low-energy constants are found we use them to give the phenomenological implications 
on the forthcoming Kn experiments. We review the experimental assumptions of [fL7| finding some 
of them only borderline compatible with the theoretical expression. We also show the form-factors 
for different values of se and present the partial- waves predicted from the full result. This provides 
insight in the relevant assumptions that the new experimental analysis should make. 
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The results were then used to obtain predictions for the tt-tt threshold parameters and the it-it 
phase-shifts. Good agreement with the existing experimental results was obtained as is visible from 
Table [| and Figs. [l4| and |l5|. We also warned the reader about the uncertainty in the relation 
between the theoretical knowledge of Kn decays and tt-tt scattering, reflected in the unknown 
0(p e ) corrections to the relation between the two- and three-flavour low-energy constants. This is 
the reason we did not include tt-tt scattering data in our fit as was done in |Bj . We do not expect 
large changes in the main results when this is done. 

The variations of the vacuum expectation values, decay constants and masses in terms of the 
lowest order parameters were shown in Figs. |l6| and [Tt] (at fixed ratio m : m s ). They show, up to 
C(p 6 ), no evidence of strong infrared singularities. 
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A Some G(p ) Contributions 

In this appendix we give the explicit formulas for the 0(p 6 ) part of the form- factors F and G of 



Ki4 split in parts as in Eq. (4.3) 



A.l LI, LI x Uj and L[x One-loop Contributions 

The F LL part of F is 

F LL = 7r4/(16^ 2 ){2/3i^ + 2/9L r 2 - 31/108L§} 

+tt4{ - 256L[L r 5 + 64L r 2 L r 5 - A8L r 3 L r 5 + 256(L r 4 ) 2 + 416L r A L r 5 - h\2L\L\ - h\2L\U % 
+56(^) 2 - 64Z^ - UL\V % - 16L^} 
+m 2 m|/(167r 2 )|8/9^ + 8/27X3} 

+mlm 2 K ^Ll (2mL\ + ML r 3 - ML\ - 32L r 5 - 128L r 6 + 8Lg) + 512(Z^) 2 - 1024L^} 
+m 2 m 2 /(16^ 2 ){-2^} 

+m 2 {L r Jl28s w Ll + 32s^L r 3 + 8s w L r g - 32UL r 2 - 16UL 3 ' + 8UL r 9 - 32u w L r 2 + 8u w L^j 

+A(m 2 )( - 404/37^ + 268/3L r 2 - 13/3LJ + 334/32^ + 3l/2L r 5 - 55/6^) 

+A(m 2 K )(- 4SL[ + 124/3LJ, - 53/3X5 + 32L^ + 25L r 5 - llLg) 

+3(m 2 ) (l6L[ + 20/3L r 2 - 73/18LJ - 10L r A + 73/6^' - 7/2i^ } 

+m|-/(167r 2 ){ - 2/3LJ - 19/9^ - 55/108^} 

+m 4 K { - ML r 2 L r 5 - WL r 3 Li + 8(L r 5 ) 2 + 8L r 5 L r 9 ] 

+m 2 f m 2 /(16vr 2 ){8^} 

+m 2 K {L r 5 ( - l28s^L\ - 32s n L r 3 - 8s^L r 9 + 32t n L r 2 + l6t n L r 3 - 8UL r g + 32u v L\ - 8u^L^j 
+A(ml) (54L r 2 + 21/2L% + 4/3LI - U/3L r 5 - 55/12Lg) 
+A(m 2 K )(- A/3L[ + 130/3L'£ + l/2L r 3 + 2AL\ + 8L 7 5 - U/2L^\ 
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+A(ml) ( - 32/3LJ + 10/3L r 2 - 49/18LJ + AL\ - 8/3L r 5 - 7/4Lg) } 
+ m y{lQn 2 ){u 2 + 1/2^-6^} 
+rr%A{m*){&L\ + 2U 2 + L r 3 + 3L r 5 } 

+B{m%, ml, s^m^OLl + 32L r 3 - 80L r 4 - 46L r 5 + 48L r 6 + 24L r 8 + 6Lgj 
+B(ml , m^>Hj - 32L r 4 + 3L 9 '} 

+B(m^ , ml , sjm^ / 3 | s?r ^ _ 200L5; - 68^ + 116L r 4 + 46L r 5 - 21L^ - 9^ - 9^^ 

+B(to^, ml, s v )^m 2 K (32/3s v L r 4 - 2s 7r L^ + 2s^g + 2s^UL r 9 + 2s v u v L r 9 ^ 

+B(ml , m^, t 7r )m^{8/3L^ - 8/3^ + 14/3^ + 3L£ - Lgj 

+B(m^, m^, t«)mlm 2 K { - 8/3L r 2 - 16/3L r 3 + 8/3L r 4 + 5L r 5 - 20/3L r e - 3/2L r 9 } 

+B(ml, m%, tv)ml{s«/2Ll + t* (8/3L r 2 + l6/3L r 3 + U/3L r 4 - 3L r 5 + 3/2^) + u v /2L\ 

+B(ml, m 2 K , t^)m A K \8/3Ll + 4/3L£ - 40/3Lg - 20/3L£ - 1/2^} 

+B(ml , m 2 K , t v )m 2 K \l/2s v L r 9 + 4t,LJ + 2t n L r 5 + t v L r 9 + l/2u n L r 9 } 

+B(ml,m 2 K ,U){ - l/2a*UL$ - 1/2*^5 - l/2t v u n L r 9 } 

+B(ml, m 2 K , u„)m* \$Q/3L r 2 + QA/3L r 4 + 20/3L r 5 - 8/3L r 9 } 

+B(ml , m 2 K , u„)mlm 2 K {ml3U 2 + 32i£ + 16L£ - 8/3^} 

m\, u 7r )4/3m^|s 7r Lg + i^Lg - 10u n L 2 - 8u^U 4 - Au^L\ + 2w 7r Lg j 
+B(ml, m 2 K , u^m A K { - l6/3L r 4 - 8/3L r 5 - 2/3Lgj 

+B(ml, m 2 K , u^)m 2 K ^2/3s^Ll + 2/3t n L r 9 + 16/3u n L r 4 + 8/3u n L r 5 + 4/3u 7r L 9 ' j 

+B(ml, m 2 K , - 2/3s 7r M 7r Lg - 2/3f 7r u 7r L 9 ' - 2/3^9} 

+B(m 2 K , m 2 K , s v )mlm 2 K ^96L[ + 2AL r 3 - 96L r 4 - 24L£ + 96L£ + 48Lgj 

+B(m 2 K , m 2 K , s^m^j - 256/3^ + 80/3^ - I6L3 + A2L r 4 + 10L£ - 4L 9 '} 

+B(m 2 K , m\, s^)m 2 K s^ - 80L[ + l6L r 2 - 20L r 3 + 28L\ - 2L r 9 } 

+B{m 2 Kl m 2 K , s^^slLl + 48s^ + 24^^ + 2s\h r 9 - \^l3s^L\ - AQ/3s,UL r 2 

—As^t^L^ + 2.?^^^ — lG/Ss^u^Li — 40/3s 7r u 7I -i2 — 4s 7r u 7r i3 + 2s 7I -u 7r Lg| 

+B(m 2 v , m 2 K , U)ml{ - 4L r 3 - 28/3L r 4 - 16/9^ + 8/3L£ + 32/3L^ + 8L r s - 8/3L^} 

+B{m 2 q , m 2 K , t n )mlm 2 K /3{2/3L r 3 - A6/3L r 4 - 22/3L r 5 + 20L r 6 + 32L r 7 + 16L r s - 2L r 9 } 

+B{m 2 v , m 2 K , U)m 2 j3{ (3L r 3 + 4L^) (s„ + «,) + 43/3t 7r L^ + l2t w L r 4 - 2t v L r 5 + 10UL r 9 

+B(m 2 , m 2 K , U)m 4 K /3{ - U/3L r 3 + 20L r 4 + 38/3L r 5 + 8L r 6 - UL r 7 - 28L£ + i^} 

+B(m 2 n , m 2 K , U)m 2 K /3{ (-3L r 3 - L 9 ') (s„ + u^j + 47/3^ + 24^L^ - 4^L£ + 2t 7T L 
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+B(m 2 ,m 2 K ,t w ){ - 4L r 4 m 2 K m 2 - 2L 7 4 m 2 T m 2 - 2L£m^} 

+B(m 2 ,m 2 K ,t 7T ){ - l/3s v t v L r 3 - s^U 9 - T j2,t\U z - 4-^9 - l/M^u^Ll - i^u^Lgj 
+B(m 2 , m 2 , s n )m^56/9L r 4 + W/3L r 5 - 16L r 6 + 64L r 7 + 24L r 8 - 2Lj} 
+B(m 2 , m%, s«)mlm 2 K { - 22A/9L r 4 - 16/3L£ + 64ig - 64L£ - L r 9 } 
+B{m 2 ,m 2 7] , s^mlml^lU/ZLl + 32/9L r 3 - 2AL r 4 - 4L r 5 } 
+B(m 2 , m 2 , s^m^j - As n L r 4 - 2s^L r 5 + s v L r 9 + t v L r 9 + u^Lgj 
+B(m 2 ,m 2 , s n ){l6s n Llrn? K + m 2 (^- 2As v L{ - As^L^j } 
+Si(mJ , , s>J { - 8L^ + 32L^ + 40L£ + 64L£ + 32L r 8 - 8Lg} 
+BiK,m^, S7T )mlm 2 K { - 32L r 2 - 8L r 3 + 64L r 4 - 4L r 9 } 

+B 1 (ml, ml, s^m^j - 368/3s„L[ + 80/3s n L 2 ~ - 32s^L r 3 + 24s n L r 4 + 8/3s w Lg + 16t n L r 2 
+4t v L r 3 + 4t v L r 9 + 16u v L r 2 + 4 UjT L r 3 + 4^} 

+BiK,m^, s n )m 2 K ^128/3s w L r 2 + 32/3s v L r 3 - 32/3s n L r 4 - 2/3s w L r g ^ 
+BiK,m', s.) {200/34^ + 12s 2 w L r 2 + 68/3slL r 3 + 2/3^^ - 64/3*^5 
-16/3s w UL r 3 + 2/3s w t v Ll - 64/3s v u w L r 2 - 16/3s w u„L r 3 + 2/3s 7r w 7r ^| 
+B 1 (ml, m 2 K , t w )m 4 w { - 8L[ - 20L r 2 - 18L r 3 - 46/3L r 4 + 16/3L£ - 4/3Lgj 
+Bi {ml, m 2 K , t n )mlm 2 K /3{l6L; + 8L r 2 - 2L r 3 + 56L r 4 + 36L r 5 - 76L r 6 - 48L r 8 - 4Lg} 
+B 1 {ml, m 2 K , t v )ml/3\8s v L\ + As^L r 2 + 5s n L r 3 + 2s v L r 9 - 2U*L\ + 12t v L r 2 + 39t*L$ 
+22UL r 4 - 20t w L r 5 + 10UL r g + 8u v L\ + Au w L r 2 + 5u v L r 3 + 2u n L^ 
+B 1 {ml,m 2 K ,t«)m i K /3\8L r 1 + AL r 2 - AL 3 - AOL r 4 - 20L r 5 + AOL r 6 + 20L r 8 — Lgj 
+B 1 (ml, m 2 K , t w )m 2 K /3{ - 8s^L\ - As v U 2 + S7T L r 3 + s n L r 9 - 8ULI + AUL r 2 
+17UL r 3 - 4UV 4 - 2UL r 5 + 5UL r 9 - 8u v L[ - 4u w L 2 + u n L r 3 + u^igj 
+B 1 (ml,m 2 K ,t w ){8s n t w L r 1 +4s w t 7r L r 2 - s„t w L r 3 - As^L r 9 - 8t\L r 2 - V3tlU 3 
— 4^Lg + 8t^u^L\ + At^u^L 2 — t v u w L 3 — 4t 7r u 7r Lg|/3 
+Bi(m2,m|-, Ux )TO^{ - AOL[ - 20V 2 - 20L r 3 - 56/3L^} 
+Bi(ml,m 2 K ,Ui V )mlm 2 K { - 32/3LJ - 16/3L r 2 - 16/3LJ - 16/3L^} 
+B 1 (ml,m 2 K ,u n )ml/3^32s n L r 1 + 16 S7I L r 2 + 16s v L r 3 + 32UL[ + l&UL r 2 + 16t n L r 3 
+48u v L[ - A8u v L^ + 2Au 7T L 1 3 + 8u^L r 4 - Au w L r 5 + 4u 7r Lg| 
+Bi(m2,m|-,u 7r )TO^/3{ - 40i^ - 20L r 2 - 20L r 3 + 24L^} 
+B 1 (ml,m 2 K ,u w )m 2 K {l6/3si T L r 1 + 8/3s«L r 2 + 8/3 S7r L r 3 + 16/3t w L{ + 8/3t w L r 2 
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+8/3t n L r 3 + l6/3u v L[ - 32/3u n L r 2 + 8/3u v L r 3 - 56/3u v L r 4 - l6/3u n L r 5 + 2/3^^} 
+Bi{m 2 n ,m 2 K ,u 7! )l/3^y — IQs^u^Vl — 8s^u v L 2 — 8-s^u^L^ — 2s 7r ii 7r L 9 — \%t^u^L\ 
-8Uu n L r 2 - 8Uu n L r 3 - 2Uu w L r 9 + 2±u\L\ + h2u\U 2 + \2u%U z - 2ulL r 9 } 
+Bi(m^, ml,U)mlm 2 K { - lO/3L r 4 + 20/37^} 
+Bi(m^,m2,^)m^{-20/3L^-10/3^ + 40/3L^ + 20/3Lg} 
+B 1 {m 2 K i m 2 K , 3^1^256/33^11 - 80/3s n L r 2 + 16s v L r 3 - 36s n L r 4 - 2s n L r 5 + 2s n L r 9 } 
+Bi(m 2 K , m 2 K , ST T )m 2 K ^32/3sT T L r l - GA/Ss^L^ + 16s^L r 4 + s^Lgj 

m 2 K , s^){ - 80/3 S lL[ - 200/3^ - 22«£L§ - s 2 n L r 9 + 16/3M*£J 
+40/3s 7r t 7r i2 + 45,^7,- — s-rrtirLg + 16 ^s^u^L^ + AQ/3s^u^L 2 + As n u w L 3 — s v u w L 9 ^ 
+Bi(m^, ro*, t^mj {l/3LJ + 1/3LJ - 8/3L£ - 8/3L r 8 } 

+B 1 (m 2 K , 171% t^[mlm 2 K (2h/3L r 4 + 7L r 5 - 20/3Lg) + m 2 n ( - 3t T LJ - 3t w Lg) } 
+B 1 {m 2 K , m 2 v t n ){m A K (46/ZLl + 2/3L^ - 8/3L£ - 4/3Lg) - 6t v L r 4 m 2 K } 
+B 1 (m 2 rj ,m 2 K , U)ml{25/3L r 3 + 55/3L r 4 + 16/9L r 5 - 8/3L r 6 - 32/3L r 7 - 8L r 8 + 17/3^} 
+B l {m 2 nl m 2 K , tv)mlm 2 K /3{ - 56/3L r 3 + 73L r 4 + 22/3L£ - 20L r 6 - 32L r 7 - l6L r s - 5/2L 
+Bi(m 2 ,m 2 K ,U)ml{ - 19/9s n L r 3 - 17/6s v L r 9 - 79/9t w L r 3 - 3t w L r 4 + 2UV 5 - 35/6UL 
-19/9^^ - 17/6^^1 

+B 1 (m 2 ,m 2 K , t„)m 4 K /3{ - 7/3L r 3 - 7AL r 4 - 38/3L^ - 8L r 6 + 64L r 7 + 28L r s - 11/2^} 
+B 1 {m 2 , m 2 K ,t^ 2B r 4 m\m\ + AL\m 2 K m 2 q + 2L\m\ } 

+B 1 {m 2 n m 2 K , t„)m 2 K {25/9s n L r 3 + 11/68*1% - 41/9t„L r 3 - 6t n L r 4 + 4t„L r 5 + l/3t«L r 9 
+25/9u n L r 3 + ll/6u w L r 9 } 

+B 1 (m 2 ,m 2 K , U){ (l/3L r 3 + 3/2Lg) (s„ + u^U + ^t\L r 3 + 3/2^} 
+Bi(ro*, m 2 v , s„){(ml + m>^) (32/3L^ + S/QLj) } 

+B 1 (m 2 ,m 2 !7 s„)ml{ - 112/3s w L^ - 32/9^ + 24s v L r 4 + As,L r 5 - 16/3^ 
-4/9UL r 3 - ml3u^U 2 - 4/9^i^} 
+B 1 (m 2 ,m 2 l , S7T ){24 S 2 7[ L r 1 + \2s\L\~ + 8s 2 n L r 3 } 
+B 21 (ml, ml,s„){ (ml + m>^) (32^ + 8U 3 ) } 

+B 21 (ml , ml , s„)[ml ( - 32s^V 2 - 32/3s„L r 4 - 40/33^1 + 8/3 S7I L r 9 - 16t^L r 2 

-U n L r 3 - \Qu^V 2 - Au^Ll} + m 2 K ( - U/3s^L\ + i/3s n L^ } 

+B 21 (ml , ml, S7r )4/3{444z- + g^LJ + 16s 2 w L r 3 - a%Ll - s n t n L r 9 - s«u„L r 9 } 
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+B 21 (ml , m 2 K , t^)ml[AQL\ + 20L r 2 + WL r 3 + IQL\ + 9L\ - 3L r 9 } 

+B 21 (ml , m 2 K , t^)mlm 2 K ^ - 112/3^ - 56/3^ - 10/3^ - 8L r 4 - 9L r 5 + 3/2L§} 

+B 2 i("4 ; m 2. ) t 7r ) m ^/3{ - 32s^L\ - 16s„L r 2 - 8s n L r 3 + s w /2L r 9 + 48t n L r 2 + 36t„L r 3 

-52t n L r 4 - 23t n L r 5 + 19/2t n L r 9 - 32u n L[ - 16u n L r 2 - 8u n L r 3 + u n /2L r 9 } 

+B 21 (ml, m 2 K , t 7r )ml{ - 8/3LJ - A/3L r 2 + A/3L r 3 - 8L r 4 + 3/2Lj} 

+B 21 (ml , m 2 K , ^Wk^/ZskLI + 16/3s v L r 2 + 2/3s r L r 3 - 3/2s n L r 9 - 16/3t T LJ 

-20/3^ - 2/3t n L r 9 + 32/3u 7I L\ + 16/3u w L r 2 + 2/3u n L r 3 - 3/2u n Ll} 

+B 21 (ml, m 2 K , U){ - 5/6s w t w L r 9 - 8t\L\ - 44/3#L£ - 34/3^ _ 5/61*1$ - 5/6^^} 

+B 21 (ml, m 2 K , u„){ (20m* + 16/3m>^) (2LJ + L 2 ' + Lj) } 

+B 21 (ml , m^, U7r ){ (l6/3ro* + 8/3™^) (s w + **)(- 2LJ - U 2 - Lj) } 

+B 2 i(m£ ; TO 2. f U7r )| m ^ (40/3L1' + 20/3^ + 20/3i^) - 8u\L\ + A/3u 2 n L r 2 - Au\U 3 } 

+B 21 (m 2 K ,m 2 K , S7r ){64/34^ - 56/3^ + 2^^} 

+B 21 (m 2 K ,m 2 v ,t 7 r){m$ (9L r 4 + 9L^ + m\m 2 K (9U 4 - 9L^) + m* (- 3^ - 3^) } 
+B 21 (m 2 K ,m*,t v ){m A K (- 18Lj) + ( - 6^) } 

+B 2 i(m^ m^, t,){mj ( - 13/3L§ - 9L^ - 3Lg) + m£m|- (56/9LJ - 9L r 4 + 3/2Lg) } 

+B 21 (m 2 ,m 2 K ,t n ){ml((lO/9s n + 4< w + 10/9^)^ - t n L r 4 - 4/3t n L r 5 

+ (3. S7r + 5t* + 3u*) L r 9 /2^j + m% (j/9L r 3 + 18L r 4 + 3/2L^ +m 2 K (- 16/9^ 

-3/2s w L r 9 - 2t n L r 4 - 8/M n L r 5 - t v L r 9 - W/9 Ull L r 3 - 3/2 U7V L^ } 

+B 21 (m 2 ,m 2 K ,t v ){ - l/2 S7T UL r 9 - 5/3i^ - 1/2^ - l/2Uu n L r 9 } 

-B 21 (m% m 2 v s„){(ml + m>^) (32/3LJ + 8/9Lj) } 

+B 21 (m 2 , m 2 v , s n )ml{l6/3t w L r 2 + 4/9UL r 3 + 16/3u*L£ + 4/9^^} 

+B 21 (m 2 v ,m 2 , S7T ){-12slL r 2 -AslL r 3 } 

+B 22 (ml, ml, S ^)ml\212l3U 2 + 56/3L r 3 - 32j3L\ - 40/3L^ + 8/3L§} 
+B 22 (ml , m^, S7r )m^{l28/3L^ + 32/3L^ - 64/32^ + 4/3i^} 

+B 22 (ml, ml , s^^UGs^Ll - 128s^L r 2 + 2A S7r L r 3 - 4s w L r 9 - 64t w L r 2 - At v L r 3 - 4UL r 9 
-64 U7r L r 2 - 28u v L r 3 - 4u^ig|/3 

+B 22 (ml, m 2 K , tn)ml/3{ - 32L1' - l\2V 2 + 88L r 3 - l96L r 4 - l04L r 5 + 32Lg} 
+B 22 (ml , m 2 K , t v )m 2 K [AL T 3 - 88/3L\ - 20/3L^ + 16/3L^} 

+B 22 (ml, m 2 K , - 16/35^^ + 40/3 S7r i 2 ' - 16/3 S7r i?; - l6/3s n L r 9 + 6A/3t n L{ 
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-2M n L r 2 - At n L r 3 - 16/3UL% + 32/3u w L[ + A0/3u w L r 2 - A/3u w L r 3 - 16/3u w L r g } 
+B 22 (ml, m 2 K , u^ml /3(l60L; + 80L r 2 + 80L r 3 + 16L^ + 16m 2 K (2L\ + L 2 + L r ^j } 
-B 22 (ml, m 2 K , U7r )8/3{ (2L\ + L 2 + (s„ + t w ) + 10u w L{ + 3u w L r 2 + hu v L r 3 } 
+B 22 {m 2 Kl m 2 K , s„){ml (32/3L\ + 368/31^ + 32L^ + m 2 K {32/3L\ + 8Q/3L r 2 + 8L r ^j } 
+B 22 {m 2 Kl m 2 Kl s^^^s^L^ - 2Q8/3s^L r 2 - 12s n L r 3 - 16/3t n L[ - 6A/3t w L r 2 - 6t w L r 3 
-16/3^^ - 16/3u n L r 2 - 2u w L r 3 } 
+B 22 (m 2 K , m 2 v , t n ){m 2 n ( - 30L r 4 - 30L^) - 60m 2 K Llj 
+B 22 (m 2 v , m 2 K ,t w )ml{268/9L r 3 + 26L r 4 - 4/3L£ + lOLg} 
+B 22 (m 2 n , m 2 K ,t n )m 2 K {8/9L r 3 + 52L r 4 - 8/3V 5 + 5L r 9 } 

+B 22 (m 2 ,m 2 K ,t 7r ){ - 22/3s n L r 3 - 5s n Ll + 2Q/3t^L r 3 - 5t v Ll - 16/3u v L r 3 - bu^V^ 
+B 22 (m 2 , m% s n ){ml (ll2/3L r 2 + 136/9^) - 24^2 - 
+B 31 (ml, ml, S7T )(ml + m^) { - 128/3^ - 32/35^} 
+B 31 (ml, ml , s v )8/3{ - 6^ - 3^ - 3s 2 n L r 3 + (8L r 2 + 2L r 3 ) (u„ + t^a*} 
+B 31 (ml, m 2 K , t*){m% ( - 32L{ - 16L$ - 4Lj) + m>| ( 32L r + mL r + 4L ^ j 
+B 31 (ml, m 2 K , t*)m* { (8L[ + AL r 2 + Lg) (s„ + u^j + 8t^L r x + At n L r 2 - t w L r 3 } 
+B 31 (ml, m 2 K , U)m 2 K { - 4^ - i w /3 + u^j (2L\ + l£) - (s„ + U/3 + u^)L$} 
+B 31 (ml, m 2 K , t„){ ( - 8L[ - 4L r 2 + L r 3 ) ( S?r + u^U/3 + (8L[ + 4L r 2 + 
+B 31 (ml,m 2 K ,u 7r )(ml +m 2 K /3^{ - lQu n L\ - 8u K L 2 - 8u v L r 3 } 
+B 3 i(ml, m 2 K , u 7r )8/3|2s 7r w 7r Li + s^u^L^ + s^u^L^ + 2i 7r M 7r Li + t n u w L 2 + Uu^L^ 
+B 31 (m 2 K , m 2 K , s v )(ml + m 2 K ^j { - 32/3s n L[ - 8Q/3s^L r 2 - 8 S7T L r 3 } 
+B 3 i(m 2 K ,m 2 K ,s )4/3{ (AL\ + lOL r 2 + 3L^j (u + u v ) s^} 
+B3i{m 2 K ,m 2 1 ,t w ){ml(^ - s v -t w - u^L r 3 + m 2 K (s v - l/3t v + u^U^ 
+B 31 (m 2 K ,m 2 v ,t w ){4L r 3 (m* - m^m^ ) + (1/35^-^ + 1/3^)^} 
+B 32 (ml, ml , s„){ml (teL\ - 160/3^ + 32/3i^) - m\ (256/3L 2 ' + 64/3i^) } 
— B 32 (ml, 771^,3^)32/3165^27^ + 7s v L 2 + As^L^ — At^L^ — t w L 3 — Au^L^ — ^£3} 
+B 32 (ml , m 2 K , t*){ml (96^ + 48L^ + 8Lg) - m^(32/3i^ + 16/3i^ + 8/3Lj) } 
+B 3 2(to^, m^, - 64/35^^ - 32/3 S7r L 2 ' - 4/3 S7T L r 3 + 112/3^^ + 56/3^5 
+10/3t 7r L?J - 112/3^^ - 56/3u w L r 2 - 10/3^^} 

-B 32 (ml, m 2 K , U7r ){l6/3(3m2 + m \ + s v +t n - u^j (2L[ + L 2 + Lj) } 
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+B 32 (m 2 K , m 2 K , s v ) (ml + m|) { - 64/3L[ - 160/3X^ - l6L r 3 } 

+B 32 (m 2 K , m 2 K , s w )| - 32 S7r X^ - 80 S7r L^ - 24s w L r 3 + 80t w Ll + 104UL r 2 + 48UL r 3 

— 16u n Li + 56?i7r i^j /3 

+B 32 (m 2 K , m 2 , t w )L r 3 { - 12m 2 + 4/3m| + 8/3^ - 14/3^ + 14/3u„} 
+3(m 2 ){ - 4/3s v L[ - 16s v L r 2 - l9/3 S7r L r 3 + 55/12s n L r 9 - 27t w L r 2 - 27/2t n L r 3 
+55/12^ - 27u n L r 2 + 55/12?^} 

+A(m 2 K ){ - be/Ss^Ll - 32/3 S7r L r 2 - 10s n L r 3 + ll/2 S7r L r 9 - U/3UL r 2 + 17/3UL r 3 
+ll/2UL r 9 - 110/3^^ + ll/2 U7T L r g } 

+A(m 2 ){ - \2 Slx L\ - s^Ll - 3t n L r 2 + ll/6UL r 3 - 3u^L r 2 + 7/4^ +U+ u^L r g } (A.l) 
For G form-factor the expression reads 
G LL = m4/(167r 2 ){ - 2L[ - W/9L r 2 - 49/108X5} 

+mt{ - 1GU 3 LI + 32L r 4 L r 5 + 56(L r 5 ) 2 64X^ - 64X5X5 I6X5X5} 
+mlm 2 K /(l6ir 2 ){8/9L r 2 + 8/27X5} 
+mlm 2 K {&4LlLl - 32(L r 5 ) 2 - 128L r 6 L r e + 8X5x5} 
+mlm 2 n /(16n 2 ){-2/3Ll} 

+ml^8s w L r 5 L r 9 + S2UL r 2 L r 5 + l6t w L r 3 L r 5 + 8UL r 5 L r g - 32u v L r 2 L r 5 + 8u v L r 5 L r 9 

+A(m 2 )( - 4L\ - 10L r 2 - l9L r 3 + 46/3X^ + 161/6X5 - 83/6X5) 

+A(m 2 K ) (32L r 2 + 17/3X5 - 8/3L r 4 + 25/3X5 - 19/3X5) 

+A(m 2 ) (8/3LI + 2/3L r 2 + 53/18X5 - 2L\ - 11/6X5 + I/2X5) } 

+m%/(16TT 2 )^2L r 1 - 7/9L r 2 - 1/108X5 } 

+^{16X5x5 + 8(X5) 2 + 8X5x5} 

+m 2 K m 2 /(16ir 2 ){8/3L r 5 } 

+m 2 K { - 8s«LlL r 9 - 32t„L r 2 L r 5 - Wt w L r 3 L r 5 - 8t v L r 5 L r 9 + 32u v L r 2 L r 5 - 8u n L r 5 L r 9 
+A(m 2 )( - 27/2L r 3 - 4/3L\ - 2/3L r 5 - 83/12X5) 
+A(m 2 K ) (aL\ + 26L r 2 + 5/6X5 + 32/3X5 - 8/3L r 5 - 19/6X5) 
+A(m 2 ) ( - 32/3X5 - 8/3X5 - 47/18X5 - 4L r A - 8/3L r 5 + l/4Lg) } 
+m^/(167r 2 ){^ + 1/6X5 -2LS} 

+m 2 {8A(m 2 )L r 1 + 2A(m 2 t )L r 2 + 3A(m 2 )L r 3 + A(m 2 t )L r 5 } 

+B(ml,m 2 K , t v )ml{ - 8/3L r 2 + 8/3L r 3 - \4j3L\ - 3L r 5 + Xg} 

+X(m 2 , m 2 K , t v )mlm 2 K {8/3L r 2 + 16/3X5 - 8/3X^ - 5L r 5 + 20/3Xg + 3/2Xg| 
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-B(m\ , m\ , t v )ml/2{s w L r g + 2t v (8/3L r 2 + l6/3L r 3 + U/3L r 4 - 3L r 5 + 3/2L^j + u v L 

+B(ml,m 2 K , t«)m A K { - 8/3L r 4 - A/3L r 5 + 40/3L r e + 20/3L r 8 + l/2Lg} 

+B(ml , m 2 K , U)m 2 K { - l/2s v L r 9 - 4t % L r 4 - 2UL r 5 - t v L r 9 - l/2u„L r 9 } 

+B(ml ,m 2 K , t 7I )^l/2s n t 7I L r 9 + 1/2^ + l/M^Lj} 

+B(ml , m 2 K , u v )m 4 n ^56/3L r 2 + QA/3L\ + 20/3L r 5 - 8/3L r 9 } 

+B(ml , m 2 K , u n )mlm 2 K {40/3L r 2 + 32L r 6 + l6L r 8 - 8/3L r 9 } 

+B(m 2 ,m 2 K , u 7r )4/3m^|s 7r ig + t^L^ - lOu v L 2 - 8u n L 4 - 4u„Xg + 2u T Lg| 

+B(ml , m 2 K , u n )m 4 K { - 16/3LJ - 8/3L r 5 - 2/3L r 9 } 

+B(ml, m 2 K , U7r )m 2 K ^2/3s^L r 9 + 2/3UL r 9 + \Q/3u^L r A + 8/3u*Ll + 4/3ii 7r ig| 
+B(ml,m 2 K ,u n )^ - 2/3s n u n Ll - 2/3t^L r 9 - 2/3ii^ig| 
+B(m 2 , m 2 K , t v )m*{4L r 3 + 28/3L r A + 16/9L r 5 - 8/3L r 6 - 32/3i^ - 8L r s + 8/3L r 9 } 
+B(m 2 , m 2 K , t v )mlm 2 K /3{ - 2/3L r 3 + AQL\ + 22/3L r 5 - 20L r 6 - 32L r 7 - 16L r s + 2L r 9 } 
+B{m 2 n , m 2 K , U)mlm 2 q \2L\ - 2L r 6 } 

-B(m 2 n , m%, t^ml/3^ (s* + u^j (3L r 3 - 4Lj) + t n ( - 43/3L^ - 12L r 4 + 2L\ - WL r 9 
+B(m 2 , m 2 K , U)m 4 K ^lA/9L r 3 - 2Q/3L r 4 - 38/9L r 5 - 8/3L r 6 + U/3L r 7 + 28/3L r s - 1/3L 
+B(m 2 ! , m 2 K , t^)[m 2 K m% (aL\ + 8L^ - 4L r 6 m*} 

+B(m 2 v , m 2 K , t v )m 2 K /3{ (3U 3 + Lj) (a„ + u^j - t„ (l7/3L r 3 + 2AL\ - AL\ + 2Lj) } 
+B(m 2 , m 2 K , t v ){l/3s w t w L r 3 + s n UL r 9 + 7 j3t\L r 3 + 4 L 9 + l/3t^L r 3 + Uu n L^ 
+B 1 (ml, m 2 K , t n )m*{8L[ + 20L r 2 + l8L r 3 + A6/3L r 4 - 16/3L r 5 + 4/3^} 
+Bi {ml, m 2 K , t n )mlm 2 K /3{ - 16L{ - 8L r 2 + 2U 3 - 56L^ - 3§L\ + 76L r 6 + A8L r s + AL 
+B 1 (ml,m 2 K ,t n )ml/3^ - 8s T iJ - 4 S7V L r 2 - 5s w L r 3 - 2s^L r 9 + 24t v L[ - 12UL r 2 
—39t w L 3 - 22ULI + 20t n L r 5 - 10t^L r 9 - 8u w L[ - Au^L r 2 - bu„L r 3 - 2^^ 
+B x (ml, m 2 K , t n )m 4 K /3{ - 8L[ - AU 2 + AU 3 + 40i^ + 20L£ - 40L^ - 20L r s + L r 9 } 
+B 1 (ml, m 2 K , t^m^/S^s^Ll + As^L r 2 - s w L r 3 - s n L r 9 + 8t v L{ - At n L r 2 - 17UL r 3 
+At 7! L\ + 2UL r 5 - HULlg + 8u^L\ + Au v L 2 - u^L 3 - M^igj 

+B 1 (ml ,m 2 K ,t v ){(s„ + U ( - 8L\ - AU 2 + L 3 + 4Lg) + t\ (8U 2 + 13U 3 + AU 9 ^ ' 
+B 1 (ml,m 2 K ,u„)ml{ - AOL\ - 20L r 2 - 2QL r 3 - 56/3L^} 
+Bi(to^ j m \ ; Ul ) m 2 m |( - 32/3^ - \<oj3L r 2 - l6/3L r 3 - 16/3^} 
+Bi(mJ, m 2 K , U7r )A/3ml { (8L{ + AU 2 + 4Lj) (a„ + t w ) 



40 



+u v (l2L[ - l2L r 2 + 6L r 3 + 2L\ -L r 5 + Lj) } 

+Bi(to^ > m ^ ; U7r ) m 4.| _ 40/3^ - 20/3L r 2 - 20/3LJ + 8Lj} 

+Si(m^, m|-, u^)8 /3m 2 K ^2s^L\ + s w L 2 + s^Ll + 2t^L\ + t^L 2 + t^L 3 

+2u w L{ - 4u n L 2 + u^L r 3 - lu v L\ - 2u^L r 5 + I/Au^Lq j 

+Bi(m 2 , m 2 K , u 7r )l/3| — Ws^u^L^ — 8s v u v L 2 — Ss^u^L^ — 2s 7r u 7r Lg 

-\U^L\ - 8Uu w L r 2 - 8t v u w L r 3 - 2t w u w L r 9 + 2<±u\L\ + h2u\L r 2 + \2u\L\~ - 2u\U 9 } 

+Bi(m 2 K ,ml, t„)\Q/z{mlm 2 K (Ll - 2^) + m A K (2L\ + L r 5 - 4L r 6 - 2L r 8 ) } 

+B 1 (m 2 K ,m 2 n t 7r )mi{ - l/3L r 4 - 1/377, + 8/3L r 6 + 8/3L r 8 } 

+Bi(m^,m^,){ro>|(- 25/3LJ - 7L^ + 20/3L[j) + 6m 2 K 6t w L r 4 } 

ro*, t*){3m% (LI + Lj) + m^/3( - 46LJ - 2LJ + 8L£ + 4L^) } 
+Bi(m* ; m 2. ) t T ) m 4 / 3 | _ 25L^ - 55LJ - 16/31,?; + 8L r 6 + 32L r 7 + 2AL r 8 - 1727g} 
+B 1 (m 2 ,m 2 K ,t 7I )mlm 2 K /3^56/3L r 3 - 13U 4 - 22/3L r 5 + 20L r 6 + 32U 7 + l6L r s + 5/2^} 
-Bi(m=,m|:,t w ){mX( 2i 4-4/3L5) + n&m*(4L; + 16/3^) + 2m%Ll\ 
+B 1 (m 2 ,m 2 K ,t n )ml{l9/9s w L r 3 + 17/6*,^ + 79/9UL r 3 + 3t*L; - 2t*L£ 
+35/6^^9 + 19/9^3 + 17/6^X9} 

+Bi(m^, m^, t^m^ / 3^7 /3L r 3 + 7AL r 4 + 38/3L r 5 + 8L r 6 - 6AL r 7 - 28L r 8 + ll/2Lg} 
+B 1 (m 2 r m 2 K ,U)m 2 K { - 25/9s n L r 3 - U/6s*L$ + Al/9UL r 3 + 6t n L r 4 - 4t w L r 5 
-l/M w L r 9 - 25/9u n L r 3 - ll/6u w L r 9 } 

-B 1 (m 2 v ,m 2 K ,t n ){ (l/3L r 3 + 3/2^) ( Sjr + u^jU + 4t 2 w L r 3 + 3/2^^} 
+B 21 (ml,m 2 K ,U)m^ - 40L[ - 20L r 2 - 10L r 3 - l6L r 4 - 9L r 5 + 3L r 9 } 
+B 2 i(ml,m 2 K ,t w )mlm 2 K {ll2/3L r 1 + 56/3L r 2 + W/3L r 3 + 8L r 4 + 9L r 5 - 3/2Lj} 
+B 21 (ml,m 2 K ,t n )ml{32/3s n L r 1 + 16/3s„L r 2 + 8/3s n L r 3 - 3/2s„L r 9 - 16t w L r 2 - 12^U 3 
+52/3ULI + 23/3t w L r 5 - 19/6t w L r g + 32/3u n L[ + \Q/3u^L r 2 + 8/3u^L r 3 - 3/2u^L r 9 ^ 
+B 2 i(ml,m 2 K , U)m 4 K {8/3L[ + A/3L r 2 - A/3L T 3 + 8L\ - 3/2L T ^ 
+B 2 i(ml,m 2 K ,t w )m 2 K { - 32/3^ - 16/3s w L5 - 2/3s n L r 3 + 3/2s n L r 9 + 16/3t v L r 4 
+20/3t 7r L^ + 2/3ULI - 32/3u w L[ - 16/3u n L r 2 - 2/3u w L r 3 + 3/2u n L r ^ 
+B 21 (ml,m 2 K ,t w ){5/2s„t v L r 9 + t\{24L\ + ML r 2 + 3AL r 3 + 5/2L^) + 5/2^^X^/3 
+B 2 i(ml,m 2 K , u„){ (20mi + 16/3m>^) (2L[ + L 2 + L r 3 ) } 

+B 21 {ml,m 2 K ,u 7 ,)lQ/3ml^ ~ 2s^L\ - s n L 2 - s„L r 3 - 2t w L[ - t v L r 2 - t v U 3 + u^Llj 
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+B 21 (ml , m 2 K , u^m^^O/SLl + 2Q/3L r 2 + 20/3^} 

-B 2 i(ml, m 2 K , u w )8/3m 2 K ^2s n L[ + s^U 2 + s v L r 3 + 2t n L[ + UL 2 + UL 3 } 
+B 21 (ml,m 2 K ,u 7[ ){ - 8u\L\ + 4/3^ - 4ujLj} 

+B 21 (m 2 K , m 2 v , U)[ - 9m% (l\ + Lj) + 9™^ TO |_ £ _ ^ + ^ + 3m 2^ ^ + ^ j 
+B 21 {m 2 K , m 2 n , i 7r ){l8m^l8L^ + 6m^L£ + m% (l3/3Lg + + 3Lg) } 
+B 21 (m 2 v m 2 K , t„)mW K { - 56/9L r 3 + 9L r 4 - 3/2L r g } 

-B 21 (m 2 v m 2 K , U)ml{ (lO/9L r 3 + 3/2Lj) + Ujr ) - ^ ( - 4L£ + + 4/3g - 5/2L 
+B 21 (m^m^i^m^ -7/9^-18^-3/2^} 

+B 21 (m 2 v m 2 K , U)m 2 K { (l6/9L r 3 + 3/2^) ( Sjr + + 2t w L r 4 + 8/3t n L r 5 + t n L^ 

+B 21 (m 2 p m 2 K , ^){l/2M.ig + 5/3^3 + 1/2*^5 + l/^u^ij} 

+B2 2 K,m^, Sw ){m* (8L£ _ 3 2L r 4 - 24L r 5 + 8Lg) + (8L^ + 4Lg) } 

+B 22 (m 2 , m 2 , s^A^As^Ll - 2s^L 2 + 2s w Ll - s v Ll - t v L r 3 - t^L\ - u^L r 3 - w^-Lg} 

+B 22 (ml,m 2 K , U)ml{ - 64/3LJ - 80/3L r 2 - 64/3L r 3 + 4/3L r 4 - 4/3L£ + 4/3^} 

+B 22 {ml , m%, t n )m 2 K \32L\ + AL r 3 - 8/3L r 4 + 20/3L r 5 + 2/3^} 

+B 22 (ml , m 2 K , t n )^16/3s v Ll + 8/3s„L r 2 + 4/3s w L r 3 - 2/3s^L r 9 + 32/3ULI + 8t„L r 2 

+16t n L$ - 2/3^ - 32/3^^ + 8/3 U7T L r 2 - 8/3u v L r 3 - 2/Su n L^ 

+B 22 (ml,m 2 K , u«)ml { 160/31^ + 80/3^ + 80/37^ + 16/3Lj} 

+B 22 {ml,m 2 K , u v )m 2 K \32L\ + 16L r 2 + 16L r 3 } 

+B 22 (ml,rn 2 K ,UT T ){ - 16/3s n L[ - 8/3s n L r 2 - 8/3s n L r 3 - 16/3t n L[ - 8/3t n L r 2 
-8/3t w L r 3 - 80/3u n L{ - 8u w L r 2 - 40/3^7^} 
+B 22 (m 2 K , m 2 K , s w )ml{32L[~ + l2L r 3 - \2L\ + 4Lg} 
+B 22 (m 2 K , m 2 K , sJrr&faLl + 12U 3 - 16L r 4 + 2L r ^ 

+B 22 (m 2 K , m 2 K , s 7r )2|2s 7r L3 — s^ig — 8t w L{ — 3t^L 3 — t n Lg — 8^^ — 3u n L 3 — u v L 
+B 22 (m 2 K , m*, t*){ml(- 6L r 4 - 6L^ - l2m 2 K L r 4 } 
+B 22 {m 2 v m 2 K , t v )ml{ - l6/9L r 3 + WL r 4 + A/3L r b + 2L^} 
+B 22 (m 2 v m 2 K , t v )m 2 K \28/9L r 3 + 20L r 4 + 8/3L r 5 + Lg} 

+B 22 (m 2 , m 2 K , - 2/3s^L r 3 - s^Llg + 16/3t w L r 3 - t„L r 9 - 8/3^X3 - u^Lg} 

+B 31 (ml,m 2 K ,U){(mi - m>^) (32^ + 16L r 2 + 4L^ } 

-B 31 (ml,m 2 K , U)ml ^4:(2Ll + L^j ^ +U+ u^j + s^L r 3 - UL 3 + u^Lg} 
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+B 31 {ml,m 2 K ,U)m 2 K {( S7t + u^j (8L[ + 4L r 2 + - 8/3t v L[ - A/3UV 2 + 1/3^} 
-B 31 (ml,ml,u v )8uv{(ml + m 2 K /3) (2L[ + L r 2 + L^j } 

+B 31 (rrii T ,m K ,Uir W{2s^L\ + s^Ll + s^Ll + 2t^L\ + t^U 2 + t^L r 3 } 
+B 31 (m 2 K ,m 2 ,t n )L r 3 ^4:( - m% + m\m 2 ^ - l/3s w t n + t 2 - l/3t w U w } 
+B 3 i(m 2 K ,m 2 v t 7r )L 3 ^ml(^s n + t n + + m 2 K (- s OT + l/3t n - u^H 
+B 32 (ml,m 2 K ,U){m 2 J64L r 1 + 32L r 2 + 12^) - m 2 K l3^L\ + 32L r 2 + 4L^j } 
+B 32 (ml,m 2 K ,U){ - 32s^L\ - Ws w L r 2 - 8s^L r 3 - \\2UL\ - 56t n L 2 
-10UL r 3 + lQu n L{ + 8u^L r 2 - 2u n L r 3 }/3 

-B 32 {ml,m 2 K ,u„){l6/3(3ml + 3m 2 K - s w -t n + u^j (2L\ + V 2 + Lj) } 

+B 32 (m 2 K , m 2 K , Sn ) (ml +m 2 K ){- 32L\ - 16L r 2 - 16L r 3 } 

+B 32 (m 2 K , m 2 K , s^^WULl + 8UL r 2 + 8UL r 3 + l&u^L\ + 8u n L r 2 + 8u^L r 3 } 

+B 32 (m 2 K , m 2 , i 7r ){8( - to 2 + m 2 K /3^L r 3 + A/3s^L r 3 + U/3UL r 3 - 2/3u*L5} 

+A(ml) (83/ Us^Ll + 27UL r 2 + 27/2t n L r 3 + 83/12t n L r 9 - 27u n L r 2 + 83/12it w i^ 

+A(m 2 K ) (l9 /6s^L r g + U/3UL r 2 - l7/3UL r 3 + l9/6UL r 9 - 110 /Zu„L$ + 19/6u„Ll\ 

+A(m 2 ) ( - l/4s„L r 9 + 27U + 3t„U 2 - U/&t n L r 3 - l/4t^L r 9 - 3u^U 2 - l/4u^L r 9 ) (A.2) 



A. 2 Cq Contribution 

Here we collect the relevant contribution coming from the tree diagram involving a 0(p e ) La- 
grangian vertex fl6[| . Each C\ is finite and as was explained in Sect. |H[ at present, for their 
numerical estimate one has to resort to models. 
For the F form-factor we obtain 

F ct = rr4j _ 3 2Cl - 32C£ - 24C 6 r - 64C? + 8C$ + 8C[ + l6C r u - 80C[ 2 - 176C[ 3 + 32C[ 4 
+8C[ 5 + 128C[ 6 - 48C[ 7 - 32C r 23 + 16C 2 r 5 + 32C£ 6 + 128Q 8 + 8C 6 r 4 + 8Q 5 + 8C 6 r 6 - 16C 6 r 7 
+32Q 8 + 16Q 9 + 16C r S3 + 32C r s4 + 8C r 90 } 

+mlm 2 K { - 16CJ + 8C r h - 40Q' - 64Cf - 24C 8 r + 16C[ + 48C[ 1 + 24C[ 2 - 96C*[ 3 - 16C[ 4 
+80CT B + 64C[ 7 - 16C 2 r 2 - 16C 2 r 3 - SC^ - 16C 2 r 6 - 32C 2 r 9 - 64C 3 r - 32C 3 r 6 + 16C r 63 + 20C r eA 
+4Q 5 - 4Q 6 - 24CJ 7 + 16Q 8 + 8Q 9 + 8Q 3 + + 12Q } 

+m^s OT {20C[ + 64CJ + 8C 3 + 8CJ + + 16C 6 r + 32Cf + 16C[ 3 + 8C 2 r 2 + 8C* 2 r 5 - 4C 6 r 4 



-4C<35 ~ 12 C66 + 8 C67 _ 32 Q8 - 16Cg 9 - 8C 83 - 16C 84 - 4C 90 | 
-m 2 



- 20C[ - MC r 2 + 8C 3 r - 4C^ - 8C r 6 - 8C r & + 8C[ 2 + 32C[ 3 + 12C 2 r 2 + 16C 23 - 20C. 
+4Q 3 - 4Q 4 - 4CSb + 6Q 6 + 16Q 7 - 16Q 8 - 10Q 9 - 12Q 3 - 16Q 4 + 8Q 8 - 10C r 90 ] 
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+mlu^ - 12C[ - 64C 2 r - 8C 3 + 36C1 - 8C[ - 16^ + 8C[ 2 + 80C[ 3 - 12C 2 r 2 + I6CJ3 
+4C£ 5 - 4Q 3 - 4Q 4 - 4Q 5 - 6CSe + 16CJ 7 - 16Q 8 - 6Q 9 - 4C& - 16Q 4 - 6Q } 
+m^{8CJ + 16C 6 r + 8C[ + 32(7^ - 8C[ 2 + 8C 2 r 2 + 16C£ 3 - 8CJ 4 + 8C 63 + 8C 64 
-4C r 66 -8C r 67 + AC r 90 ] 

+m 2 K s^Cl + 8C 3 r + 16CJ + 32Q + 32C 7 ' + 16C£ + 8C[ 2 + 32C[ 3 - 16C* 2 r 3 + 8C* 25 - 8Q 3 
-8Cg 4 + 8Cg 7 - 8C 68 - 4Cg 9 - 4C 90 j 

+m 2 K t n { - AC[ + 8C 3 + 4CJ - 8C£ - 16C 6 r - 24C[ 2 - 64C[ 3 - 4C£ 2 - 8C r 23 - 12C 6 r 3 
-8Q 4 + 6CSe + 12Q 7 + 2Q 9 + 4C& + 4C& - 10Q } 

+m^ r t* w {4C[ - 8C 3 - 4CJ - 8C[ - 32^ + 16C[ 2 + 32C[ 3 - 4C 22 - 8C£ 3 - 4Q 3 - 8Q 4 

+6Cg 6 + 12Cg 7 — 2Cg 9 — 4C 83 — 2C 90 | 

+4{ - 8C{ - 32CJ + 8CI + 4Q 6 + 8C& + 4Q 9 } 

+s^{4C[ + 32CJ - 4CI - 2Cgg - 4Q 7 + 8C& + 6Q 9 - 2Q 8 + 2C 9 r } 

8CJ - + 32C r 2 - 8C r 3 - 20C r 4 + 6C r 66 - 4C 6 r 7 + 8C r 6S + 2C r 69 - 2C r 8S + 2C 9 r } 
+tl{lC{ - ACl - 2C r 66 - 2C\ 7 - 2C r 69 - 6C r ss + 6C r 90 ] 
+t„u n { - 4CT - 8C 3 r + ACl - 6Q 6 - 12Q 7 + 2Q 9 - 2Q 8 + 2C r 90 ) 

+ul{8C r 3 -2C r 67 }, (A.3) 
while for the G form-factor the contribution is 

G ct = { - 24C 4 r - 8C 6 r - 8C 8 ' + 8C r 10 + WC^ - 16C[ 3 + 8C[ 5 + 16C[ 7 + 8C 2 r 2 - 8C£ 5 
+8Cg 3 + 8Cg 4 + 8Cg 5 - 4Cg 6 + 12Cg 9 + 8C 83 + 12C 90 | 

+mlm 2 K { - 16CI - 8C r b - 24Q - 8C r & + 16C[ + 486*^ - 24C[ 2 - 64C[ 3 + 16C[ 4 
+16CT B - 16C 2 r 2 + 16C r 25 + WC r 26 - 32C r 29 + 16Q 3 + 20C 64 + 4C r 65 
+4Qg + 8Q 9 + 8Q 3 + 12C 9 r } 

+^%{ - 4CT + 8C 3 r + 8C£ - 16C[ 2 - I6CT3 - 8Q3 - 4Q 4 - 4Q 5 + 4Q 6 
-8Cg 9 + 4C 88 - 12C 90 | 

+mlu[AC{ + 8C r 3 + YlCl + 8C r 6 + 8C r 8 - 8C[ 2 + 4C 2 r 2 + 16C 2 r 3 + 4C r 25 - 4Q 3 - 4Q 4 
-4Q 5 - 2Q 6 - 4Q 7 - 10Q 9 - 4C£ 3 + 4C 8 r 8 - 10C 90 } 

+mlu^ - 4C[ - 8CI + 20C r 4 - 8C[ - 16(7^ + 8C[ 2 + 16C[ 3 - 12C* 2 r 2 - 16C 2 r 3 
+4C 2 r 5 - 4Q 3 - 4Q 4 - 4Q 5 + 2Q 6 + 4Q 7 - 6Q 9 - 4Q 3 - 6Q } 

+m^{ - 8Q' - 8C r b - 16C 6 r + 8C[ + 32C r u - 40C[ 2 - 64C[ 3 - 8C 3 r 4 + 8C 6 r 4 + 4C 6 r 9 + 8Q 3 } 
+m 2 K s^[ - AC[ + 8C r 3 + 8C[ 2 - 32C[ 3 + 8C r 22 - 8C r 25 - 8C r 6i - AC r 69 - 8C r S3 + AC r S8 - 4CJ } 

44 



+m 2 K u{AC{ + 8C r 3 + 4CJ + 8C r 5 + 16CJ + 8C[ 2 + 4C 2 r 2 + 8C 2 r 3 - 4C 6 r 3 - 8C 6 r 4 - 2C r 66 



"6Cgg — 4Cg 3 — 2Cg | 



+m 2 K u 



ic U'Tt^ — AC\ — 8C 3 -(- 12C 4 — SCJq — 32CJi ~l~ 16C^ 2 -4- 32C]^ 3 — 4C 22 — 8C 23 

~ 8 C23 ~ 4^63 _ 8 C64 _ 2^66 ~ 2Cg 9 ~ 4 Cg 3 - 2C 90 j + - 4Cg 8 + 4Cg | 

+s 7r t 7r {4C[ - 4CJ - 2Q 6 + 6C& - 2C& + 2C£ } 
+s 7T u 7r -^ — 8C3 — 4C*4 — 2Cg 6 + 2Cgg — 2Cg 8 + 2Cg j- 
+4{ - 4C[ + 4CI + 2C& + 2CJ 7 + 2C r e9 2C r 88 + 2C r 90 } 

+t^{±C{ - 8C 3 r - 12CJ + 2Q 6 + 2C£ 9 - 2C 8 r 8 + 2C 90 } + u\ {8C r 3 - 2C£ 7 }. (A.4) 
A. 2.1 Resonance Contribution 

The aim of this section is to give an estimate, based on Resonance Saturation as discussed in Sect. 
|, of the previous combinations of 0(p 6 ) constants. The main drawback of this method is that 
it does not specify the scale at which it should be apply. A rough estimate of this uncertainty is 
obtained varying the scale \i, let's say between 0.5 and 1 GeV, and comparing the results with the 
scale at the p mass. 

Using the equations of motion for integrating out the massive fields we get for the F form-factor 



F RS - p2 



{ c li( _ 8r7 4 + 12mlrn 2 K - 4m|-) + c m c d (A:m% - 4mlm 2 K + 8mjs r - im 2 t n ^ 
+c^( - 8m* a* + 2mlu + 2m 2 K U + 4-4 - 2#) } 

+ p2 a^a r4 \fvfx(~ 4 "4 + m l(-6 m K + 2s tt +t n + 3u n ) - 2m% + m%(2s 7r + 3t n + u n )) 
b ^J2M v i. \ ' 

+gvf x (l2m* + l0m 2 m 2 K — irri^s^ + I0m 2 t n — 6m 2 u T + 2m A K — Ajt? k s v — 2m 2 K uSj 
+gvav(2m^ — m 2 7I (m 2 K + s„ — t w + u^) - m 4 K + m 2 K (s 7r + 2t„ + u^) + s^u^ — 2t^{s^ + 
+V2g v fv ( - 2mlu - rn 2 K t^ + (s n U + s^it^ + 'it 2 + t 7r ii 7r )/2^ 
+V2g v (5m 2 t n + m^t^ - s^t^ - s^u^ - Uu^j + V2f 2 \32m^ + IQm^rrv 2 ^ 
+V2f x av(^Qm 4 + 8m 2 m 2 K — S,m 2 s w — 12m 2 — im^U^ + Am^t^ - Am 2 K u^ 
W2/2(32mt + 16m>^)} 

' lH {fAl [ A ] [-ml-m 2 K +u^+ f Al {2) (iml - m 2 K - 2s„ + t») } , (A.5) 



F 2 M\ 



and for the G form-factor 



^{ 4c m(™^ m l: ~ m K) - 4c m c d ml(ml - m 2 K - - 2t^<?J i m\ + m 2 K - t^j | 
8/7?,^ + lAm 7I m K — 2m 7r t 7r — om w u w + 2m K + Am K s 7T — 2m K u 7r ) 
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+g v ay ( — 2m^. + to 2 (to|- + Ss^ + t v + u v ) + m A K - m 2 K {s^ + u v ) — 2s„-i 7r — s^u^j 
+V2g v fv ( - ml(s v + t w ) - m 2 K s^ + s 2 + (s n t n + s^u* + t\ + t n u n )/2j 
+V2gy ((2s^ + t^iml + m 2 K ) - s n t n - s^u^ - Uu^j 

+V2f x a v (Sml + 8m\m 2 K - 4m 2 r (t 7r + ti T ) + 8m 4 K - Am 2 K (2s v +t„ + u T )J 
+48V2f 2 mlm 2 K } 

+ WKp{ fA1 ^ (~ m l~ m K + «*) + Ia1 { a ] (ml + m 2 K -u)}. (A.6) 



As can be seen from Eqs. ( A.5), (|A.q ) the a xial-v ector contribution vanishes at si — and thus 
has no influence on the values quoted in Eq. ( 6. 12] ). 



B Scattering Lengths 



Using the definitions of Eq. ( |6.2l|) is rather straightforward to evaluate the threshold parameters 
a\ and b\ not displayed in 0. 
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807r 3 Fi 2520 



«,2- 1 
2 240tt 3 Flml 



277 
~840 
11 



1849 1583tt 2 1 _ 124 — 179 T 773 T T 17 T 

61 09 O3 64 — Ok H Ok 

17010 22680 84 315 378 210 3 

~193tt 2 24218 1 T 337 — 157 T 87 — - — 

6i 6 2 H 63 H 64 - 3 o 5 + 5 b 6 

7560 2835 7 420 126 14 



5607r 3 F*m2 [168 

47 

:!o2807r 3 i^m^ 1 15 



4111 37?r 2 



1134 840 
549221 41tt 2 



1 17 151- 653 — — — 

Ol H O2 63 O4 + O5 + Ok 

14 84 126 126 



5400 



18 



15 5 



(B.l) 



where 



16tt 2 F2 : 



2.3.4 



167T 01,2,3,4, 65,6 = (167T ) 6 5 , 6 



(B.2) 



For the rest of the notation and the expression for the 0.; coefficients, we refer to App. D in 0. 



C Vacuum Expectation Values at G(p 6 ) 

Here we collect the relevant formulas for the vacuum condensate in the limit m u — rrid- As can 
be seen they can be cast in terms of the simplest one-loop integral A defined in App. 
For the non-strange current we find 

(0|W|0)( 4 ^ -L{3/2A(ml ) + A(m 2 K ) + A(m 2 )/6 + 4m 2 (4Lg + 2L r s + + 32m 2 K L r 6 } , 



(0\qq\0) 



1 r 1 



Ft l (16tt) 



m 4 (47/3888tt 2 + 47/648) - to 2 to^(113/19447t 2 + 113/324) 



+m 2 m 2 (19/4327T 2 + 19/72) + m|-(19/4867r 2 + 19/81) - m 2 K m 2 (13/216tt 2 + 13/36) 



+to*(5/216tt 2 + 5/36) 



+21/8(A(m 2 )J + 7/2^1(m 2 )A(m 2 f ) + 7/12^(m 2 )A(m 2 ) 

+A(ml)(ml(-24L r 4 - \2U b + U2L r 6 + 68L r 8 + 10ff 2 r ) + H2m 2 K L r ^j + ((A(m 2 K ^ 
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-A{m 2 K )A{ml)/Q + A(m 2 ) 2 /72 

+A{m 2 K )(ml {8L r 5 + 32L r 6 + 8L r 8 + 4ff 2 r ) + m 2 K {-32L\ - WL r 5 + 128L r e + 327^)) 
+A(m 2 n )(ml/3{20/3L r 5 - 16L r 6 + 64L r 7 + 12L r 8 - 2H r 2 ) 
-m 2 K /3(32/3L r 5 - 112L r 6 + 64L r 7 ) - 8m 2 L^ 

+64mt((LJ + Ll){ALl + 2L\ + flj) 4(7^) 2 6L r 6 L r a L r & H r 2 2(L r 8 ) 2 L r 8 H r 2 ) 
+128m 2 m 2 K {{Ll - Lg)(8Lg + 2U 8 + H r 2 ) + 2L\L^ 
+m%(l024L r i L r 6 + 256L r 5 L r 6 - 1024(7^) 2 - 512L r 6 Ll) } 

+^{^ (l2C[ 9 + 20C r 20 + l2C r 21 - Q 4 ) + mlm 2 K (±8C r 21 + + m% (l6C£ + A8C r 21 ) } , 

(C.l) 

where (0|gg|0) stands for (0|mt|0) either (0| cZcZ|0) . The results for the strange quark reads 
<0|^|0)( 4 >= ^{A(m 2 K ) + A(m 2 v )/3 + m 2 (87^ - AV 8 2H 2 ) + 4m 2 K (4L r 6 + 2I7 8 + H r 2 )} 
(0\ss\0) (6) = ^4 { J^-y frV 6 + 1) ( - 7/162,74 + 13/81m>| + l/36m 2 m 2 + 4/81m 4. 
-2/3m 2 K m 2 + 17/36m*) 

+£A(ml)A(m 2 K ) + 4/3A(m 2 )3(m 2 ) + 2(A(m 2 K )) 2 + 2/9(A(m 2 )) 2 

+A(ml) (m 2 (-247^ + 88L r 6 - 20L r 8 - 1077 2 r ) + 16m 2 K (AL r 6 + 2L r 8 + 77 2 ')) 

+A(m 2 K ) (4m 2 (AL\ + l2L r 6 - 2L r 8 - H r 2 ) + 8m 2 K {-AL r A - 47^ + 20L' 6 + lOL r 8 + 77 2 r )) 

+A{m 2 ) (ml i '3(80/37^ + 8L r 6 - 1287,; - 60L r 8 + 2H r 2 ) 

+m 2 K /3(~128/3L r 5 + 160ig + 1287,; + 1447^ + 8H r 2 ) ~ 8m 2 L 4 ) 

+m A J{Ll + LI- Lg)(256Lg - 128Lg - 6477J) + 128(7^) 2 + 64L r s H^) 

+64m>| (l6L^ + *L\L\ + 2U h U 8 + L r b H' 2 - 16(7^) 2 ) 

+64m%(4(Ll + Ll){ALl + 2U 8 + H r 2 ) - AL\(AL% + AL r 8 + H r 2 ) - 4(7^) 2 - 2U 8 H r 2 ) } 
+ i{4mj (l2CT 9 + AC r 20 + 12C r n + Q 4 ) - Umlml (3C[ 9 + C r 20 - 3C r 21 ) 
+192m 4 K (c r 19 + C r 20 + C r 21 )}. 

(C.2) 

D Loop Integrals 

In this Appendix we collect for completeness some familiar formulas for the one- loop integrals. 
The two-loop integrals that appear we discussed in [^) and for the others we used the expressions 
of [^7j . We have not discussed these integrals in more detail here since we do not present any 
formulas involving them. 

Through the calculation one has to use one-loop integrals of one, two and three point functions. 
The latter disappears after mass renormalization and the use of some recursion relation. All in all 
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i J (27r) d q 2 - m 2 ' 
d d q 



i J (2ir) d (q 2 — m\){{q — p) 2 — m 2 ,) 
d d q % 



i J {2v) d (q 2 -m 2 ){{q-p) 2 -m 2 ) 
P^Bi{m\,ml,p 2 ) , 

d d q q^q v 



i J (2Tr) d (q 2 — m\){{q — p) 2 — m 2 ) 
p fJ ,PuB 21 (ml,m 2 ,,p 2 ) + g^B 2 2(m 2 ,ml,p 2 ) , 
d d q quisle 



we only have to deal with the following set of functions - in the remainder we use d = 4 — 2e. 

u 2, 1 f d d q 

A { m i) 

TDt 2 2 2\ 

B{m 1 ,m 2 ,p ) 
BJm 2 ,m 2 ,p 2 ) 



D / 2 2 2 



R ( 2 2 2\ _ 

= P^PvP a Bz\{m\,m\,p 2 ) + {p^g va + p^g^ a + p a g^v)B Z 2(jn\,m\,p 2 ) . 

(D.l) 

An expansion in e leads to the following series 

2 

A{m\) = ■^X + A{m 2 ) + eA\m 2 ) + ... , 
B l0 {m 2 ,m 2 ,p 2 ) = j^pole^ + B^imlmlp 2 ) + eB^^mj, m 2 2 ,p 2 ) + ... , (D.2) 

with A, Bij defining finite quantities and where "pole^" denotes the singular part of each of the 
Bij functions, 

pole = A , pole x = -^, pole 21 = ^ , pole 22 = ^(m 2 + m\ - ^r) , 



2 r Z1 3 r " 4 v 1 3 
4 ' r ~- JZ 24 



P° le 3i = ~r ' P° le 32 = ^r( 2m i + 4m 2 ~ P 2 ) . ( D - 3 ) 



with 

A = - + ln(47r) + 1 - 7 . (D.4) 
e 



After some simpler algebraic manipulation, the functions defined in Eq. (D.l) can be related 
to the basic integrals A{m\) and Bi(m1,m 2 ,p 2 ) through the identities 

B zl {m\,m 2 2 ,p 2 ) = -^(A(m\) - (m\ - m\ - p 2 )B 2 i{m 2 1 , m 2 2l p 2 ) - 4B 32 (to 2 , m 2 ,p 2 )) , 
B 3 2(m 2 ,mlp 2 ) = ^^-^A(m 2 1 ) + ^A(m 2 )-(m 2 -m 2 -p 2 )B22(m 2 ,m 2 ,p 2 )y 
B 2 i(ml,m 2 ,,p 2 ) = ^(^A{ml) +m 2 l B{m\,ml 1 p 2 ) - dB 2 2{m\ 1 ml,p 2 )^ , 
B22{m 2 l7 m 2 ,p 2 ) = — - - (A{m 2 ) + 2m 2 B(m 2 , m 2 ,p 2 ) 

-(p 2 +m\ - ml)B x (m\,ml,p 2 )^ , 
B 1 {mlm 2 ,p 2 ) = ^(A{m 2 )- A(m 2 ) + (m 2 ~m 2 +p 2 )B(m 2 ,mlp 2 ] 

B(m 2 ,m 2 ,0) = { ±^A{m 2 ). (D.5) 
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Notice that the inclusion of the previous identities is only done at the final (numerical) level 
in order to avoid cancellations between different terms occurring in the form-factors. Then the 
one-loop contribution is reduced to 



A(m?) 

15 f 2 2 2\ 

B(m 1 ,m 2 ,p ) 



+ 



16tt 2 

1 to 2 ln(m 2 ) — m 2 li^m 2 .) 
167T 2 m 2 — m 2 

1 



(32^ 2 ) 



A S \ 771? f , (p 2 + v) 2 

n + -r hi — o o In ■ 



A 5 



(p2 _ v yi _ A 2 



(D.6) 



with A 



n\ — m 2 , S = rn 2 + ot| and v 2 = [p 2 — (mi + tt^) 2 ]!?? 2 — (mi — 1TI2) 2 }. Similarly 
combining Eqs. (D.2) and ( |D.5| ) one can obtain the expressions of the e terms in the expansion as 
functions of A and B . A straightforward calculation leads to 

16^(m 2 ) = rn 2 [^ + I + ^! + Iln 2 (m 2 )-C71n(m 2 )] 



16n 2 7?(mlmhp 2 ) = ^ - \ + g + (C m 2 , P 2 ) + \ 



with C = ln(47r) + 1 — 7 and m 2 = (1 — x)m 2 + xm 2 — x(l — x)p 2 . 



dx ln 2 (m 2 ) , 



(D.7) 
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